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1.0  INTRODUCTION 


Several  topics  relating  to  conformal  arrays  have  received  detailed 
attention  in  the  last  few  years.  These  topics  include  pattern  synthesis, 
types  of  appropriate  radiating  elements,  and  their  impedance  and  mutual 
coupling  characteristics,  and  wide  angle  scanning  of  linear  arrays  located 

on  planar  and  conical  surfaces.  

Synthesis  techniques  for  the  far  field  patterns  from  current  distri- 
butions on  planar  surfaces  are  well  known.  Since  there  are  no  proven 
techniques  for  synthesizing  antenna  patterns  from  arrays  on  conical  sur- 
faces, two  different  approaches  have  been  pursued: 

1.  An  equivalence  principle  was  used  to  determine  the  distri- 
bution of  sources  on  a cone  to  produce  a prescribed  pattern. 

Using  this  approach,  planar  arrays  are  replaced  by  sources 
on  a conformal  surface.  Thus,  it  is  possible  to  generate 
patterns  using  well  known  synthesis  techniques  for  planar 
surfaces.  The  equivalent  sources  are  continuous  functions 
on  the  cone.  The  practical  validity  of  the  method  depends  on 
the  accuracy  with  which  the  continuous  functions  can  be 
approximated  by  physically  realizable  radiators.  Equivalent 
source  distributions  were  computed  at  discrete  points  on  a 
conical  surface  to  approximate  the  patterns  of  a specific 
planar  array. 

2.  The  computation  of  approximate  patterns  from  discrete 
radiators  judiciously  positioned  on  the  cone  was  carried  out. 

An  optimization  was  performed  to  obtain  element  weightings 
that  yield  the  highest  signal-to-noise  ratio  in  the  beam 

pr  inting  direction  neglecting  the  effects  of  mutual  coupling. 

Uniform  weighting  schemes  were  also  used  for  beams  in  the 
region  about  the  cone  axis.  This  second  approach  has  been 
termed  a he  .ristic  approach. 
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In  the  Final  Report  on  Contract  000 1 9-72 -C -02 1 2,[  1 ] the  results  from 
the  above  two  approaches  were  reported.  In  the  equivalence  principle 
approach,  excitations  and  radiation  patterns  were  computed  for  two 
particular  beam  pointing  directions.  These  patterns  were  then  compared 
to  the  far-field  pattern  that  would  have  been  produced  by  the  generating 
planar  array. 

In  the  use  of  the  heuristic  approach  several  layouts  of  radiating 
elements  were  considered.  Parameters  such  as  variations  of  sidelobe 
levels,  grating  lobes  and  beamwidths  were  investigated. 

In  conjunction  with  the  use  of  the  above  two  approaches,  computer 
programs  were  written  for  the  approximate  computation  of  the  fields  and 
patterns  of  different  configurations  of  elements  on  the  surface  on  the  cone. 
The  program  used  in  the  heuristic  approach  was  appended  in  the  first 
Quar^drly  Report  [2], 

In  the  present  report,  'he  computation  of  an  exact  element  pattern 
is  < -scribed.  The  report  begins  with  the  general  expressions  for  two 
potential  functions  representing  the  modal  fields  of  arbitrary  apertures 
on  a conducting  cone.  Expressions  for  the  electric  field  components  are 
then  derived  from  the  potentia'  functions  for  circumferential  and  radial 
slots  with  particular  slot  excitations.  The  field  expressions  are  the  basis 
for  a computer  program  which  computes  radiation  patterns  from  both  cir- 
cumferential and  radial  slots.  The  appropriate  approximations  for  the 
L egendre  functions  and  the  spherical  Bessel  functions  used  for  their  com- 
putations are  summarized  in  S^etrtm~-3  of  this  report. 

The  computer  program  is  a modification  and  extension  of  one 
described  by  Pridmore - Brown  t^T^nd  has  been  executed  on  the  IBM  370 
computer  system,  now  operational  at  the  Hughes  Aircraft  Company.  The 
program  computes  the  special  functions  in  DOUBLE  PRECISION,  thus 
giving  the  required  accuracy  of  the  computations.  The  complete  program 
together  with  a brief  description  of  its  use  is  incli^fed  in  Appendix  A.  A 
second  program,  useful  in  studying  the  exact  element  modal  patterns  at 


minimum  cost,  is  included  in  Appendix  B. 
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Severa.  test  cases  were  excuted  on  the  program  in  order  to 
check  its  correctness.  In  each  test  case,  the  results  compared  favor- 
ably with  identical  cases  of  a circumferential  slot  considered  by 
Pridmore  - Brown  and  Stewart  [4], 

A comprehensive  analysis  of  the  modal  series  patterns  is  pre- 
sented in  Section  4 of  this  report  for  a slot  39  radians  from  the  tip. 

The  computed  patterns  are  then  compared  with  previously  measured 
patterns  [5],  Aiso  pres>  nted  are  patterns  of  circumferential  and  radial 
slots  at  30  radians  from  the  cone  tip  for  comparison  with  the  previous 
slot  location. 

Suggestions  for  further  investigations.  Section  5,  conclude  the 

report. 
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2.  0 RADIATION  FIELDS  OF  CIRCUMFERENTIAL  AND 
RADIAL  SLOTS  ON  A CONE 

Expressions  for  the  radiation  fields  from  Blot-excited  conical 
antennas  were  first  obtained  by  Bailin  and  Silver.  ^ The  fields  were 
expressed  in  terms  of  normal-mode  expansions  using  the  orthogonal 
properties  of  the  eigenfunctions.  Although  some  corrections  appeared 

r *7 1 

in  the  literature1  modifying  the  first  results,  the  expressions  for  the 

fields  need  additional  examination.  More  recently  Pridmore -Brown 

r 4l 

and  Stewart1  1 presented  a more  rigorous  calculation  of  the  O-polarized 
electric  field  component  of  circumferential  slots.  Their  numerical  cal- 
culation is  based  on  expressions  derived  from  integral  transform  methods. 

In  this  section,  complete  expressions  of  two  potential  functions 
representing  the  modal  fields  will  be  presented  for  both  the  circumfer- 
ential and  the  radial  slot.  The  electric  field  components  will  then  be 
obtained  from  the  potential  functions  and  will  be  available  for  calculation 
of  radiation  patterns  of  the  individual  slots  or  arrays  of  slots  of  any  de- 
sired configuration. 

Figure  1 shows  typical  circumferential  and  radial  slots  and  the 
conical  geometry.  The  usual  spherical  coordinate  system  centered  at 
the  cone  tip  is  associated  with  the  structure.  The  cone  axis  coincides 
with  the  Z-axis  of  the  rectangular  coordinates  (x,  y,  z)  associated  with 
(r,  0,  d).  Primes  denote  coordinates  of  the  apertures  on  the  structure. 

The  two  scalar  potential  functions  in  question  are  n and  n , 
TE  and  TM  to  the  radial  direction,  respectively.  For  an  assumed  time 
dependence  of  exp  jjutj  both  satisfy  the  wave  equation: 

(V2  + k2)  n = 0.  (1) 

The  electric  and  magnetic  field  components  are  then  conveniently  given 
by  superposition  of  TE  and  TM  partial  fields  obtained  from  the  potential 
functions  as  follows: 
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Figure  1.  Conical  Geometi  { and  Slot  Coordinates 
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For  an  arbitrarily  shaped  source  excitation  the  expression  for  the  n 
is  given  by  Bailin  and  Silver.  (There  is  (-)  sign  in  their  expression.  ) 
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where  r>,  r<  symbolizes  the  larger  ana  smaller  of  the  coordinates  r,  r' 
respectively.  P^(cos  0)  is  the  associated  Legendre  function  and 
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are  the  spherical  Bessel  and  spherical  Hankel  functions,  respectively. 

Jj,  is  the  Bessel  function  of  the  first  kind  and  order  v and  is  the 

Hankel  function  of  the  second  kind  a d order  v . 6 is  the  Kronecker 

om 

delta  function 
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whil"  the  V.  are  the  roots  of  he  equation 
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The  function  fj(r\  4')  in  equation  (4)  is  the  electric  field  source  excitation 
in  the  r direction.  Following  the  general  procedure  described  by  Ballin 
and  Silver  where  the  bound  ry  conditions  are  applied  in  conjunction  with 
the  Lorentz  reciprocity  thiorem  to  the  modal  fields,  an  expression  for 
nTE  is  obtained.  The  TE  scalar  function  with  arbitrary  source  excita- 
tion is  given  by, 
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and  V is  the  free  space  wave  impedance  { t ' Vi  art  roots 
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and  f.,  (r',  6')  is  the  E , component  of  the  source  excitation.  All  other 
quantities  are  as  defined  for  the  TM  case.  It  is  interesting  to  observe 
that,  unlike  the  TM  mode,  the  TE  mode  is  excited  by  both  the  r-  and  d- 
components  of  the  source  excitation  in  the  aperture.  Conversely,  the  r- 
component  of  the  source  excitation  excites  both  the  TM  and  the  TE  modes. 
It  should  also  be  noted  that  the  expression  in  Equation  (9)  differs  from  the 
expression  given  by  Bailin  and  Silver  [7]  in  the  expression  for  r (r,  r') 
and  a common  factor  of  1/sin  0q. 

The  modal  functions  and  the  corresponding  field  components  from 
particular  aperture  configurations  on  the  cone  structure  with  appropriate 
field  excitations  are  given  in  Sections  2.  1 and  2.2  below. 


2.  1 Circumferential  Slot 

A narrow  circumferential  slot  of  width  2w<<Ais  considered  as 

shown  in  Figure  1.  The  narrow  slot  has  an  azimuthal  length  d = 2 d a sin  0 

o o 

where  2 d is  the  azimuthal  angle  subtended  by  the  slot,  0 is  the  cone 
o ° 7 o 

exterior  angle  measured  from  the  axis,  _a  is  the  distance  from  the  cone  tip 
to  the  center  of  the  slot  and  \ is  the  operating  wavelength.  The  slot  is 
assumed  to  be  excited  by  a voltage  Vo  across  its  center  and  has  a resulting 
field  across  it  in  the  r direction  given  by 
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Substituting  the  expression  of  E^r',  4')  in  Equations  (4)  and  (9)  for 
i (r1,  4')  and  integrating  over  the  primed  coordinates,  we  obtain  the 
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set  equal  to  zero  in  the  TE  case,  as  only  the  Er  component  of  the  source 

excitation  is  assumed.  Thus, 
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where  f.(r,  a) 


(15) 

r < a 


Hie  raduition  fields  are  readily  found  by  substituting  the  expression  of  D 
and  n corresponding  to  r>  a into  Equations  (2)  and  (3).  On  using  the 
asymptotic  expansion  for  the  spherical  Hankel  functions. 
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the  following  expressions  are  obtained  for  the  electric  field 


components : 
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It  is  pointed  out  that  the  O-dependence  is  represented  by  the 

Legendre  function  P(>m  j ^ (cos  ®)  instead  of  P™  (cos  0)  as  originally 

derived.  This  choice  simplifies  the  evaluation  of  the  Legendre  function 

the  digital  computer  and  is  preferred.  Also,  except  for  a constant 

t.ictor  of  1 / tt  , the  expr es s ion  f or  E agrees  with  that  reported  by  Pridmore- 
[4l  W 

Lrown  and  Stewart1  1 derived  by  the  Kontrorovitch -Lebedev  transform 

method. 


2.  2 Radial  Slot 


A narrow  slot  of  width  2 w«circumference  is  assumed  to  be  posi- 
tioned along  a generatrix  of  the  cone  and  extends  from  r^  to  r^  as  shown  in 
Figure  1.  It  is  further  assumed  that  the  length  of  the  slot,  (r^  - r^),  is 
such  as  to  allow  definition  of  the  ends  of  the  slot  by  constant  d' . The  slot 
is  excited  by  a voltage  resulting  in  an  electric  field  in  the  (^-direction 
given  by 


E 


6 


v’  g(r')S(d') 
o r ' s in  0 


o 

The  function  g(r')  describes  the  variation  of  the  source  excitation  in 

the  r’  direction.  Clearly,  since  fj  ( r ' , d')  = 0,  i.  e.  , no  source  excitation 

in  the  r' -direction,  the  TM  mode  of  the  potential  function  is  not  supported. 

The  TE  mode,  on  the  other  hand,  is  obtained  only  from  the  contribution  of 

f7(r',  d')  in  Equation  (9).  Substitution  of  Equation  (24)  into  Equation  (9)  re- 

TE 

suits  in  the  following  expression  for  II  for  the  radial  slot: 


TE 
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(r) 


co  oo 


• 0 (2 l'  + 1 ) Pyi  (cos  0)  cos  m d 

17  ~i  dZPrn(cos  0) 

™-°  ‘-1  Sine  r(l  + a l 
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f ^K{ 

Jr , 1 


om  did  Q 
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(kr'<)h  , (kr>)dr' 
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0 = Q 
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The  electric  field  components  are  obtained  from  Equation  (2) 
operating  on  Equation  (25)  giving: 


' -j(kr-V4)  oo 

-E>VF 

m=  1 


m sin  m 4 


1 


r sin  0 


sin  0 0 


‘4 


V ' e ’j(k  r ~K/4)  °°  I 

- E 

m=0 


cos 


m 4 


i sin  0 


(1  + & ) 4 

om 


where 

P„ 


oo 

E 

i=l 


CD 

E .i 

i=  1 


P^m- 1/2  (cos  0) 
'iJ  ^2p;mi/2(cos  0) 


J , (k  r') 


d " d 0 


d p-m,/7(cos  0) 


r 2 

• / &iil)  1a. 

0 = 0O  r!  v k r' 

t 

V = V. 


. a 0 % : - i /2 

. 1'.  . . 1 
J 1 


a 2 p;Ti  j2  <c°9  °) 
a ► a 0 


/ 


0=0  1 
o 


V = V . 

1 


2 J.(kr') 

iiiDli 

r'  yfT~ 


and  the  v 1 are  the  roots  of  the  equation 
-m 

a P , /,  (cos  0 ) 
y -1 /2  o 
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3.0  COMPUTATIONAL  PROCEDURES  FOR 
SPECIAL  FUNCTIONS 


A brief  description  of  the  procedures  used  in  the  computation 
of  80me  functions  appearing  in  the  far  field  expressions  is  given  below. 


i.  1 


Computation  of  Legendre  Functions 


From  the  expressions  of  the  field  components  given  in  Section  2, 
n that  the  calculation  of  radiation  patterns  requires  the  evaluation 
tV  function  P^vcos  0)  and  its  derivatives  with  respect  to  r and  0.  As 
ted  previously,  we  have  chosen  to  work  with  the  function  P*m  (cos  0) 
her  than  PJ,m(cos  0).  The  two  functions  are  related  byL 


Pvm(cos0)  cos(mr)  = 


r (v-t-  m + 1 


r (v 


m 


+ 1) 


p:m(cos  o) 


(31) 


where  T (k)  is  Euler's  gamma  function.  Using  the  above  relationship,  it 
is  easy  to  show  that  the  ratio 

— P™  (cos  0) 

gQ  * 

— — P^*‘  (cos  0 ) 
e v v o 


e.in  be  replaced  b 


d 


8 
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p;m(cos  o) 


— Pjm(cos  0 ) 
d V v ' o 


in  the  expression  of  the  functions  n™  and  H • From  equation  (31)  it  is 
1 so  evident  that  the  roots  of  P™  (cos  Pq)  are  also  roots  of  P,'™  (cos  PQ). 
The  same  holds  true  for  the  derivatives  of  the  two  functions  with  respect 
to  0.  Derivatives  of  the  Legendre  Functions  with  respect  to  P are  com- 

r o i 

puted  by  the  recurrence  formula1 
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m 


(32) 


»>  = 


cos  0 T^m,  , „m+l, 

~--n-  q-  Pw  (cos  0)  + Pv  (cos  0) 


d m f 1 Ol 

derivatives  Py  are  computed  numerically  using  the  formula1  1 

•dV£,'-i)  ■ ik  {«-0»  - 8f(Vl>  t 8f(vz)  - 

, „r4,  y = v + i h,.  h is  a selected  increment  in  V,  and  l(v.)  are  the  values 

v‘  1 ° m 1 

0l  the  Legendre  functions  of  order  m and  degree  Py  (cos  0). 

Thus  far  we  have  given  expressions  for  computing  the  various 

derivatives  of  the  Legendre  functions  in  terms  of  the  functions  themselves. 

The  computations  of  these  functions,  however,  require  extreme  care,  as 

valid  approximations  must  be  used  for  the  different  ranges  of  the  variables 

‘ 9 

0 and  v . For  relatively  small  values  of  0 such  that  Q < rr:  where  0 is 

fill 

expressed  in  radians,  the  formula1 


P;m(cos 


m oo 

0,  _ / 1 - cos  0 \Z~  Y"  / nn  T (t-t  n t 1 ) 

\ 1 + cos  9 ) r ( ^ - n + 1 ) T (m+n  + 1 ] 

' ' n = o 


.1  1 n,n 

h - 2cos  0) 


n ! 


(33) 


is  suitable  for  computations.  As  before,  I'  (k)  is  the  usual  Euler's  gamma 

function.  The  computation  of  P™  (cos  0)  is  then  obtained  from  Equation  (31). 

However,  for  v large,  (v  > 100),  the  ratio  ■*'  m t— U - is  replaced  by  its 
2m  1 v p — rn  » l ) 

limiting  value  v 

9 f 12l 

For  large  values  of  0 such  that  0 > YZ+Z — we  use  t^ie  exPans’onl  J 


P;m(cosO)  - rl'-m*  !)£  (-if 

k = o 


(1  /2  ■)  m)k  (1  /2  - m)k 
k!  (2  sin  0)k  (•'  + 3/2). 


(34) 


sin  [(i/+  k + 1 /2)  0 + n (1  /2  k + 1 /2m  t 1 /4)] 


where  (b)^  - The  above  expansion  gives  satisfactory  results  for 


all  values  of  * and 


T> 


< 0 < 5 


For  0 outside  this  range,  the  series  is 


asymptotic  and  its  accuracy  depends  on  the  value  of  V 

15 


This  last  expression 


s uscd  for  numerical  computations  in  the  neighborhood  of  V=  70  + 2m  and 
t|ic.  desired  functions  are  obtained  through  backward  recursion. 


y 2 Computation  of  Spherical  Bessel  Functions 

The  computation  of  the  fields  involves  the  spherical  Bessel  function 

(2) 

j (x)  and  the  spherical  Hankel  function  hj,  (x).  Since  these  functions  are 
related  to  corresponding  cylindrical  functions  by  Equations  (5)  and  (6),  the 
cylindrical  functions  are  computed  instead.  For  the  Hankel  function, 
the  asymptotic  expansion.  Equation  (lb),  has  been  used. 

In  the  accompanying  computer  program.  Appendix  A,  the  Bessel 
functions  are  computed  in  DOUBLE  PRECISION  using  computer  system 
library  routines.  If  such  library  routines  are  not  available  or  accessible, 
high  speed  computational  subroutines  for  these  functions  may  be  written 
from  appropriate  approximations.  These  approximations  are  given  in 
Appendix  C. 

The  computation  of  the  radiation  field  from  a slot  along  a generator 
of  the  conical  surface  requires  the  ev  luation  of  the  integral 


J„  (kr) 


d r 


A series  solution  of  the  indicated  integral  when  g(r')  is  a polynomial  in 
r'  is  obtained  from  the  formula  [13] 


4. 

/ J„  (t)  dt 
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* + M -f  1 
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GO 

y (*-+  2k  + 1) 

k -o 


r + k) 

n-— ?--3+  k) 


J<'+2k+l  (z) 


with  Real  part  (4  t 1)>0.  Alternatively,  the  indicated  integral  may 
be  evaluated  with  sufficient  accuracy  by  known  numerical  integration  tech- 
niques. The  latter  method,  using  Simpson's  "3/8  Rule",  has  been  chosen 
for  the  present  program.  Either  method,  however,  requires  the  evaluation 
of  the  Bessel  function  (x). 


lb 


4.0  ANALYSIS  OF  COMPUTED  ELEMENT  PATTERNS 
FOR  SLOTS  ON  A CONE 


The  convergence  of  the  modal  aeries  representing  the  far  fields  of 
the  circumferential  and  radial  slots  was  examined  by  looking  at  the  pattern 
change  in  the  sequential  summation  of  higher  order  modes.  In  each  case, 
patterns  were  calculated  for  every  10°  in  azimuth  and  every  5°  in  elevation 
for  both  the  0 -polarization  and  the  d-polarization.  However,  elevation  cuts 
were  plotted  only  at  d = 0°  and  d = 1 80  degrees.  These  two  cuts  have  been 
plotted  on  an  extended  6 scale  to  give  a single  view  of  the  elevation  cut. 

Since  the  cone  under  study  has  a 10°  half  angle,  the  position  of  0 = 80°  was 
selected  for  azimuthal  cuts.  This  corresponds  to  the  broadside  of  the  conical 
surface.  Figure  1 shows  the  relative  position  of  slots  in  the  coordinate 
system. 

The  patterns  to  be  shown  subsequently  have  been  normalized  to  the 
largest  value  computed  for  both  polarizations  for  the  particular  modal  sum 
or  individual  mode.  For  convenience,  we  shall  refer  to  the  individual  mode 
patterns  by  the  small  letter  "mM  and  to  the  sum  patterns  by  the  capital  letter 
"M". 

4.  1 RADIATION  PATTERNS  OF  CIRCUMFERENTIAL  SLOTS 

A slot  of  a zimuthal  length  d = X /Z  located  atka  = 39  radians(k  = 2iT/A.) 

from  the  tip  of  the  10°  half  angle  cone  was  selected  for  pattern  computations. 

This  particular  location  was  selected  because  it  corresponds  to  a case  of  a 

circumferential  slot  at  8.  15  GHz  and  9 inches  from  the  cone  tip  previously 

r 1 7 1 

investigated  experimentally.1  1 The  computed  results  will  be  compared 
with  the  available  experimental  results  where  appropriate.  The  slot  is 
assumed  to  be  excited  by  a voltage  across  its  center  and  has  a resulting 
field  across  it  in  the  r direction  given  by  Equation  ( 1 2 ).  Near  the  tip  of  the 
cone  this  excitation  function  may  be  represented  by  a few  terms  of  a cosine 
series  with  appropriate  coefficients,  A greater  number  of  terms  are  re- 
quired for  adequate  representation  of  the  same  function  when  the  slot  is 
located  many  wavelengths  away  from  the  tip.  It  will  be  seen  later  that  there 


• s a correlation  between  the  number  of  expansion  terms  in  the  excitation 
function  and  the  number  of  terms  required  for  convergence  of  the  modal 
series  representing  the  fields.  However,  this  correlation  is  not  easily 
determined  beforehand.  The  dominant  TE^  mode  of  a rectangular  wave- 
guide may  be  used  to  excite  the  slot. 

Figures  2 and  3 show  the  O-polarized  pattern  of  the  lowest  two 
modes,  m = 0 and  m = 1,  respectively,  for  4 - 0°  and  4 - 180°.  As 
expected,  the  modal  patterns  for  4 = 180°  are  identical  to  the  4 = 0°  patterns, 
since  the  variation  of  the  fields  in  the  ^-direction  is  purely  sinusoidal.  The 
pattern  corresponding  to  the  m = 0 mode  shows  significant  variation  at  broad- 
side with  extrema  at  9 = 0°  and  9 = 9^  = 170  degrees.  The  peaks  at  9 = 12° 
and  9 = 30°  of  this  particular  mode  are  particularly  noteworthy.  In  contrast 
to  the  m - 0 mode  pattern,  the  m = 1 mode  pattern  of  Figure  3 exhibits  a 
relatively  small  and  smooth  variation  at  the  broadside  region.  The  maxi- 
mum now  occurs  at  9 = 0°  with  a pronounced  first  minimum  at  9 = 15° 
followed  by  a local  maximum  at  9 = 25  degrees.  Also,  the  pattern  is  about 
1 dB  down  at  9 = 9q. 

The  total  or  sum  9-polarized  pattern  of  the  first  two  modes,  M = 1, 
is  also  shown  in  Figure  2 for  4 = 0°  and  for  4 - 180°  together  with  the  M = 2 
pattern.  The  addition  of  the  second  mode,  m = 1,  to  the  first  mode,  in 
effect  dominates  the  resulting  pattern.  The  effect  of  the  m = 0 mode,  how- 
ever, is  still  evident  in  the  region  near  the  cone  axis  with  a -3.  5 dB  value 
at  9 = 1 degree.  The  forward-back  difference  varies  from  1 dB  at  9 = 1°  to 

9 dB  at  9 = 80°  and  to  14  dB  at  9 = 9 . 

o 

Higher  order  mode  sum  patterns,  M = 3 through  M = 5,  for  both 

4 = 0°  and  4 = 180°  are  shown  in  Figure  4.  In  each  case,  the  sequential 

summation  of  a higher  order  mode  in  the  sum  series  effectively  decreases 

o * 

the  variation  at  the  broadside  region  with  the  average  level  of  the  4 - 0 
patterns  at  -9  dB  for  M = 2 and  -4  dB  for  M = 4.  Addition  of  the  m = 5 
mode  raises  the  pattern  level  only  slightly  at  broadside  while  in  the  vicinity 
of  9 = 40°  and  9 = 150°  changes  of  about  -1.  5 dB  are  noted.  The  next 
significant  changes  in  the  pattern  are  seen  in  the  M = 7 modal  sum  shown 
in  Figure  5 for  4 - 0°.  At  the  same  time,  the  4 = 180°  pattern  shows  that 
the  pattern  level  decreases  with  increased  modes  at  the  broadside  region 
while  maintaining  the  level  of  about  -5.  5 dB  and  -15.  5 dB  at  9 = 1°  and 
0 = 0 , respectively. 
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The  addition  of  the  next  five  modes,  m = 8 through  m = 12,  to 
the  M = 7 total  pattern  gives  the  total  pattern  of  M = 12.  The  patterns 
for  this  case  are  also  shown  in  Figure  5.  Comparison  of  this  pattern 
with  the  M = 7 total  pattern  shows  that  the  addition  of  these  five  modes 
results  in  a decrease  of  the  pattern  level  of  less  than  .4  dB  in  the  range 
of  0 = 40°  to  0 = 140  degrees.  Similarly,  the  4 - 180°  pattern  for  this 
case  shows  only  a slight  decrease  in  the  same  range.  There  is  no  signi- 
ficant change  to  either  pattern  outside  this  range. 

Azimuthal  patterns  for  the  broadside  angle  0 = 80°  were  plotted 
for  both  0-  and  flf-polarizations.  The  *e  are  shown  in  Figures  6 through  8 
for  the  0-polarization  and  in  Figure  9 for  the  ^-polarization.  From 
Figures  6 and  7 it  is  seen  that  the  0-polarized  pattern  changes  considerably 
with  the  addition  of  higher  order  mode  almost  at  all  azimuth  angles.  The 
M = 7,  0-polarized  pattern,  however.  Figure  8,  is  down  -5.  8 dB  at  4 - 0° 
and  decreases  almost  monotonically  to  about  -26  dB  at  4 - 100°  and  remains 
below  that  level  for  the  remaining  azimuthal  region.  The  M = 7,  ^-polarized 
pattern  is  below  -23  dB  for  all  azimuth  angles,  Figure  9.  Smoother  patterns 
are  observed  in  both  polarizations  of  the  M = 12  case  as  compared  to  corres- 
ponding pattern  of  the  M = 7 case. 

To  study  further  the  effects  of  the  higher  order  modes  to  the  total 
pattern,  the  m = 13  mode  was  added  to  M = 12  pattern.  The  resulting  M = 

13  sum  pattern  in  elevation,  is  shown  in  Figure  10.  The  azimuthal  patterns 
for  the  M = 13  case  have  been  included  in  Figures  11  and  12.  It  is  seen 
that  there  are  no  noticeable  changes  in  either  the  elevation  or  the  azimuthal 
patterns  with  the  addition  of  the  m = 13  mode.  It  should  be  pointed  out,  how- 
ever, that  the  normalized  0-polarized  pattern  of  the  m = 13  mode,  has  its 
maximum  at  0 = 90  degrees  rather  than  at  broadside.  This  is  a consequence 
of  the  0-dependence  of  the  associated  Legendre  functions  and  their  deriva- 
tives. This  also  explains  why  only  the  broadside  range  of  the  pattern  is 
changed  with  increased  higher  order  modes. 

The  measured  E-plane  pattern  of  the  circumferential  slot  fed  from 
a half -height  X-band  waveguide  is  shown  in  Figure  13  together  with  the  M = 

12  computed  pattern.  The  two  patterns  agree  very  well  near  the  cone  tip 
and  for  the  broadside  region  where  both  patterns  are  quite  uniform.  The 
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Figure  6.  9-Polarized  Total  Patterns  of  A /2  Circumferential  Slot 
for  0 80°,  M = 0,  1,  and  2 
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measured  pattern  drop  off  to  the  rear  is  a consequence  of  shadowing  at 
the  base  of  the  cone  by  the  absorbant  material  in  which  the  cone  was  set. 
Thus,  it  is  concluded  that  about  13  terms  of  the  modal  series  are  required 
to  correctly  represent  the  radiation  from  this  slot  configuration.  Further- 
more, the  tip  diffraction  effects  are  generally  accounted  by  the  lower  order 
modes  of  the  series. 

Elevation  patterns  were  also  computed  for  a similar  slot  located 
at  ka  = 30  radians  from  the  tip  of  the  cone.  These  three  patterns  corres- 
ponding to  M = 4,  M = 5 and  M = 6 are  shown  in  Figure  14.  It  is  pointed 
out  that  the  two  cuts  shown  here  are  for  4 - 0°  and  <4  = 90°.  Although  addi- 
tional modes  are  required  to  describe  the  radiation  characteristics  of  this 
case  adequately,  the  figures  do  shown  the  effect  of  tip  diffraction  as  the  slot 
is  moved  along  a generatrix  of  the  cone.  The  first  minimum  of  the  4 - 0° 
patterns  now  occurs  at  0 = 25°  from  the  tip  of  the  cone  as  compared  to  0 = 
20°  for  the  ka  = 39  position  (Figure  5). 

4.  2 RADIATION  PATTERNS  OF  A RADIAL  SLOT 

In  the  present  section  we  consider  the  narrow  slot  along  a generatrix 
of  the  cone  which  was  described  in  Section  2.2.  For  the  present  study  the 
length  of  the  slot  has  been  chosen  as  one -half  wavelength  at  the  operating 
frequency.  The  cone  characteristics  and  the  operating  frequency,  as  well 
as  the  location  of  the  center  of  the  slot  remain  the  same  as  for  the  circum- 
ferential slot  described  in  Section  4.  1.  The  slot  is  excited  by  a voltage  Vq 
resulting  in  an  electric  field  in  the  ^-direction  given  by 
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o r'  sin  0 ' ’ 


It  is  interesting  to  follow  the  convergence  of  the  modal  series 
of  the  radial  slot  in  the  same  manner  that  was  done  for  the  circumfer- 
ential slot.  We  first  look  at  the  m = 0 and  m = 1 mode  patterns  of  the 
d-polarized  components  at  4 - 0°  depicted  in  Figures  15  and  16.  The 
effect  of  the  cone's  tip  on  the  modal  patterns  is  clearly  seen  by  the  be- 
havior of  the  patterns  near  the  axis  of  the  cone.  This  effect  continues 
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Figure  15.  «i-Polarized  Total  Patterns  of  \ /Z  Radial  Slot  for 
4>  - 0°,  oS  = 180°,  M = 0,  1,  and  2 
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to  persist  also  in  the  total  patterns  of  M = 1 through  M = 5,  shown  in 
Figures  15  and  17.  It  is  evident  from  the  Figures  that  the  total  patterns 
are  dominated  by  the  m = 1 mode.  Both  the  patterns  for  4 = 0°  and  4 = 
180°  exhibit  a local  maximum  at  about  20  degrees  from  the  axis.  This 
maximum  moves  toward  the  broadside  region  with  increased  number  of 
terms  of  the  modal  series.  Also,  the  higher  order  total  patterns  are 
increasingly  narrower.  The  higher  order  patterns  of  M = 6,  7,  and  12 
are  shown  in  Figure  18. 

The  M = 12  for  4 = 0°  pattern  shows  a minimum  of  -18  dB  at  0 = 1° 
and  has  its  maximum  at  exactly  broadside.  The  pattern  for  4 = 180°  of  the 


same  case  dips  to  a minimum  near  0=15  and  then  increases  to  about  -14 


d • at  0 = 100  degrees.  Both  patterns  decrease  rapidly  near  0 = 0q.  The 


irregularities  noticed  in  all  the  previous  total  patterns  are  not  present  at 
all  in  the  M = 12  patterns. 

Azimuth  patterns  for  the  ^-polarized  component  and  the  0-polarized 
component  for  the  radial  slot  are  shown  by  Figures  19  through  23.  As  it 
is  seen  from  theFigures,  all  of  the  low  order  sum  patterns  for  both  polariza- 


tions exhibit  considerable  variation  in  azimuth.  This  can  be  viewed  as  an 


indication  of  the  slow  convergence  of  the  modal  series  and  the  requirement 
for  additional  terms  in  the  series  representation.  On  the  other  hand,  the 
M = 12,  d-polarized  total  pattern  shows  a smooth  transition  from  0 dB  at 


4 - 0°  to  -13  dB  at  4 - 140°  while  maximum  of  the  9-polarized  pattern  occurs 


near  broadside  at  about  -17  dB. 


From  the  above  considerations,  it  appears  that  more  than  eight  but 
less  than  fourteen  terms  of  the  modal  series  are  required  to  represent  the 
fields  of  this  slot  configuration  correctly.  The  upper  limit  was  further 
checked  by  looking  at  the  M = 13  patterns  and  comparing  them  with  the  M = 
12  patterns.  These  patterns  are  shown  in  Figures  24  through  26.  There 
were  no  noticeable  changes  in  either  the  elevation  or  the  azimuthal  patterns. 
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Figure  18.  «J-Polarized  Total  Patterns  of  \ /Z  Radial  Slot  for 
4 = 0°.  4=  180°,  M = 6,  7,  and  12 


ilb:  Basiss 
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Figure  25.  ^-Polari 


'I 

Elevation  patterns  of  the  same  slot  located  at  ka  = 30  radians  from 
the  tip  of  the  cone  were  also  computed  and  are  shown  in  Figure  27.  As 
in  the  circumferential  6lot  for  this  position,  the  <i  - 0°  and  the  d = 90°  cuts 
are  shown.  The  maximum  of  the  M = 6 pattern  occurs  at  0 = 72°  instead 
of  the  expected  broadside,  0 = 80°.  This  is  because  of  the  insufficient 
number  of  modes  in  the  truncation  of  the  modal  series.  All  other  aspects 
of  the  pattern  are  the  same  as  for  the  ka  = 39  position,  indicating  that  tip 
diffraction  has  little  effect  on  the  radiation  pattern  of  the  radial  slot.  A 
greater  number  of  terms  must  be  taken  in  the  modal  series,  however,  as 
the  slot  is  moved  farther  away  from  the  tip  of  the  cone.  Also,  these  higher 
order  modes  affect  mainly  the  broadside  region  of  the  element  pattern. 

Although  we  do  not  have  experimental  data  for  the  slot  described 
in  this  section  for  comparison  at  this  time,  the  above  patterns  do  look  very 
much  like  the  patterns  of  a slot  on  a cylinder  tipped  10  degrees  from  its 
axis.  In  further  studies  of  the  problem,  the  above  results  will  be  verified 
experimentally. 


-Polarized  Total  Patterns  of  \/2  Radial  Slot  for 


5.0  RECR  M MEN  DA  T IONS  FOR  FURTHER 
INVES'l  CATIONS 


The  analysis  of  the  computed  element  patterns  from  circumfer- 
ential and  radial  slots  on  a sharply  tipped  cone  presented  in  Section  4 of 
this  report  is  only  a partial  one.  For  more  complete  understanding  of 
the  ra  liation  characteristics  of  such  el  merits  additional  investigative 
worl  must  be  performed.  Specifically,  the  developed  computer  program 
tor  the  exact  element  patterns  may  now  be  used  to  compute  patterns  of 
slots  near  the  tip  of  the  cone  and  very  far  away  from  the  tip.  The  results 
, n then  be  analyzed  so  as  to  separate  the  contribution  of  tip  diffraction 
from  the  total  pattern.  Such  a procedure  not  only  aids  in  understanding 
the  radiation  characteristics  but  also  will  result  in  more  economical 
pattern  calculation  as  discussed  in  the  next  paragraph. 

The  present  analysis  has  shown  that  more  modes  are  required  in 
the  truncation  of  the  modal  series  for  element  locations  far  away  from  the 
tip  of  the  cone.  The  computation  of  a higher  order  mode,  however,  requires 
considerable  more  computer  ti  le,  for  the  index  of  the  v series  in  the 
field  expressions  is  a function  of  both  the  order  of  the  mode  and  the  location 
of  the  element  relative  to  the  tip  of  the  cone.  If  the  contribution  of  tip 
diffraction  is  separated  from  the  total  pattern,  many  repetitive  computations 
would  be  eliminated  and  conical  array  problems  will  be  solved  more  econom- 
ically. Such  an  approach  will  allow  the  c oi  putation  of  patterns  to  verify  the 
equivalence  principle  pattern  synthesis  technique!^ 

The  computed  element  patterns  of  the  circumferential  slot  relatively 
close  to  the  tip  of  the  cone  have  shown  good  agreement  with  measured  patterns. 
The  experimental  verification,  however,  needs  to  be  accomplished  for  other 
element  positions  and  orientati ' 'ns . In  addition,  the  present  analysis  technique 
could  be  extended  so  as  to  take  into  consideration  mutual  coupling  effects. 

The  results  could  then  be  compared  with  similar  results  obtained  from  the 
geometrical  theory  of  diffraction. 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  CONICAL 
ARRAY  PATTERN  ANALYSIS 


A FORTRAN  IV  version  of  the  program  used  to  compute  exact 
element  patterns  of  circumferential  and  radial  slots  on  a cone  is  pre- 
sented at  the  end  of  this  appendix.  Specifically,  the  program  uses  the 
IBM  370  computer  system  presently  available  at  the  Hughes  Aircraft 
Company  facilities.  The  field  expressions  of  Sections  2.2  and  2.  3 are 
the  basis  for  the  program.  The  two  different  types  of  slots  are  con- 
sidered concurrently  when  they  are  located  at  the  same  position  on  the 
conical  surface  or  individually  when  located  at  different  positions. 

The  input  data  to  the  program  in  a specified  format  specifies  the 
type  of  slots  to  be  considered,  their  location  on  the  conical  surface,  and 
other  parameters  necessary  for  the  computation  of  the  field  patterns. 

The  program  computes  both  modal  and  total  radiation  patterns.  The 
modal  patterns  provide  the  necessary  information  for  properly  truncating 
the  modal  series,  while  the  total  pattern  information  is  stored  on  tape  for 
subsequent  use  in  array  analysis.  The  input  parameters  required  by  the 
READ  statement  are  defined  as  follows: 


THETAO 


KA 


THETAI 

THETAF 


0q,  supplement  of  the  half -angle  of  the 
cone  in  degrees . 

position  of  center  of  element  from  the  tip 
in  radians. 

initial  value  of  spherical  angle  0 in  degrees, 
final  value  of  spherical  angle  0 in  degrees. 
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DTH 


increment  of  the  spherical  angle  0 in  degrees, 
final  value  of  azimuthal  angle  6 in  degrees. 


DPMI 


increment  of  the  azimuth  angle  6 in  degrees. 


Program  always  evaluat'  s pattern  for  6 = 0. 
first  value  of  the  modal  series.  Normally  M - 0. 
If  M 4 0,  the  lower  modes  are  not  included  in  the 


computations . 


MCR 


switch  specifying  the  type  of  slot  to  be  con- 


sidered. MCR  = 1 specifies  both  circum- 
ferential and  radial  slots  individually  at  same 
position.  MCR  = 2 specifies  only  circumferential 
slot.  MCR  = 3 specifies  only  radial  slot. 


MVM 


upper  limit  of  modal  series  for  circumferential 


MVR 


upper  limit  of  modal  series  for  radial  slot. 


Some  of  the  remaining  primary  parameters  of  the  program  are 


defined  as  follows: 


counter  and  identifier  of  different  oases  con- 


sidered. 


number  of  terms  in  the  v - series  required 


for  convergence. 


PHIN 


array  of  values  of  the  azimuthal  angle  6' 


THETA 


the  spherical  angle  0. 


array  of  values  of  the  spherical  angle  0. 


maximum  number  of  devalues. 


NPHI 


maximum  number  of  d-values. 


A 


LX  total  number  of  points  of  field  evaluation. 

LM  DO  loop  index.  LM  = 1 outputs  patterns  of 

circumferential  slot.  LM  2,  outputs  pattern 
of  radial  slot. 

'additional  relationships  between  the  program  variables,  printed 
.put  and  the  field  component  equations  are  presented  for  convenience  in 
hie  A - 1 . 

In  addition  to  printed  output,  the  program  also  provides  a record  of 
.11  computed  field  component  values  on  magnetic  tape  for  further  utilization 
In  a separate  program.  The  desired  length  in  inches,  the  labeling  informa- 
tion and  the  initial  and  increment  values  for  both  axes  must  be  specified  in 
the  argument  list  of  the  CPLOT3(.  . . . ) plotting  routine.  Azimuthal  cut  plots 
.ire  obtained  only  if  JS  / 0 for  0 = ST(JS)  through  0 = ST(JL)  degrees  in  steps 
..f  JSTEP.  The  program  terminates  with  the  completion  of  the  computations 
r < • r the  last  data  card. 


CONICAL  AR-fAV  . 

I)  I Mr  NK  ! ONI  NT,  Of  THF  B J , R E S , A R , ft  R , N U ARRAYS  MUST  BE  GREATER  PR 
, j 1 1 A I !()  I N D I X "MX",  THE  NUMBER  UF  RnPTS  OF  l.E  OF  NDRF  FUNCTION. 
COMPLEX  E C T ( 8 , 665  ) , ECF  ( R,665),ERT(  8,665)fEPF(  8,665) 

COMPLEX  ACT186),ACF(85),ART(85),ARF(P5),AR(120,3),BRI120,3) 

COMPLEX  ECTW, ECFM, EKTM, ERTM 

complex  xrN,xNP,Psn,PS,PsnR,PSP,CFi,rF2,zr.R0,r)j 

DOUR  LF  PRECISION  8 J , R E S , B I , NO  I 99 , 2 ) , NO  1 , NU2  , XNU,  XNUPL,  P,  PD,  P I ,0 

CT'IPLF  PPECISIOM  SNT  ,CST  ,P  AP.PTR  , P IP2,  T TO,  PI 

01  Mf  NS  I ON  8 J (99)  ,PFS  (9r<  ,2  ) , PHI  N I 20  ) ,CMP<  20)  , SMP  ( 20 ) , MSI  15) 

01  MENS  ION  F ( 2,  15)  ,KAS(  15)  , SI665)  ,SF  (665)  , SMT  (20)  ,SMF(  20)  ,ST  ( <,0 ) 

01  MF  NS  I ON  ST1(40),ST2(40),ST3(40),ST4(40),SP(40),SFP(40) 
EQUIVALENCE  ( AR,S ) , I BR, SF) 

EQUIVALENCE  ( B J , S M T ) , ( 8 J ( 2 1 ) , S Mf  ) 

REAL  KAfKAL,KAS 

DATA  PI, ZERO, OJ/3. 1415926535 89 79 ,(0.,0.),(0.,1.)/ 

DATA  JS,JL,JSTFP/17,  17,1/ 

MCR  IS  A PARAMETER  LOR  TYPE  OF  SLOT.  MCR  = 1 TREATS  CIPClJMF. 

AND  RADIAL  CASF  AT  SAME  POSITION,  KA  DISTANCE  IPC)M  TIP. 

^C  R r 2 TREATS  CIRC.  ONLY,  WHILE  MCR  = 3 TREATS  RADIAL  ONLY. 

DATA  JS, JL , JSTEP/ . . ./  SELECTS  THE  THETA  -CUI  PLOTS. 

FORMAT  ( 1TTI  , 49H  T HE  T AO  KA  THFTAl  THTTAF  OTH  "HI  OPHI  dk, 
l 7X , I 4HM  MCR  MVM  MVR  ) 

HO  FORMAT  (2X,2F7.3,5F6.2,F8.5»4I4) 

111  FPRMAT(3E15.5) 

112  FORMAT (2X,4HMCR=,  12) 

113  FORMAT (2X,11HMCR,LM,KCR=,315> 

115  FORMAT ( 1H1, 'CONE  ANGLE  =*,F6.1,'  DEGREES,  KA  =',F8.3,’,  M =',I2) 

1 16  FORMAT ( 1 HQ • 5 HMOOF  I 1//3X , 1 HN,8X,2HNU, 19X , 1 IHCP  Fr?CIENT,20X, 
15H0PDNU, 12X, 12 H BESSEL  FUNC T ) 

117  FORMAT  ( 1 H 1 , 5HMUDE  , I l , 5X , 4HMCR= , I 2 ) 

M8  FORMAT (2X,2( 4H  NO  ( , I 3 , 1 H , , I 2 , 2 H ) = , F 1 4 . 0 ) ) 

1 2 U FORMAT ( 1H0,  • THETA  = • , F 6 . 1 , ' DEGRF F S •/ / 3 X , 1HN , 4X , 1 4HQF -L E GE N TONCT, 
1 1 1 X , 17H THETA-PARTIAL  SUM  , 10X , 1HN.4X, 14H LEGENDRE  FUNCT, 

1 1 1 X , 1 7HPHI  - PART  I AL  SOM) 

121  FORM  at ( 2 ( I 4 , 4X  , F 1 3 . 5, 5X , 1m ( , E 1 3. 5 , 3H  , ,E13.5,IH))) 

130  FORMAT ( I4,F14.9,E20.8,E16.8,2E20.8,F 16.8) 

140  FORMAT { 1H1 , 30X ,*R ADI  AT  ION  PATTERN, • ,F4.0, 'DEGREE  CONE'// 

12X, 1 1 HE  XC I T A T I ON , 2 5HE  L FM.  TYPE(CIRC=l,PAO=2»*.I3,5X, 

1 IPHE-CCMPI TH= l ,PH=2)  = , I 3/ ( 3(4X ,4HKA  =,F7.2,2X,2HM=, I 3, 2X, 

1 3 FC  = , F 7.4)  ) ) 

141  FORMAT  (1H+-,50X,,KA  = • , F 7 .2  , ?X  , • M = • , I 3,  ?X  , • FC  =',F7.4/ 

IIF24.2,! 8,F13.4,F17.2, I8,F13.4) ) 

1 A 3 FORMAT ( // 1 X , 10HPHI ,TH( TA=,F5.1,16F7.1/19X,17F7.1)) 

170  F0RMATIE6.1  ,3X  ,17F7.1/(9X,17F7.1  )) 

171  FORMAT  1 2 X , 3 ( 1X,3HAR( , 13, 1H,, I2,2H)  = ,2E13.5)  / ) 

172  FORMAT ( 2 X , 3HM Q = , I 1 ,2X,E 13.5,5X,1H( ,E 13.5.3H  , , E 1 3 . 5 , 1 H ) , 
l 3X  , 2E  1 o . 5 ) 
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. j 3 fOR,*AT(?X,4E15.5,4X,2HS<,I3,?H)  = ,F10.3 ,2X  ,2  HOB  ) 

! U)  F '1R«M!  1HO,  ' TIU  T A =•  ,Ffe.  I , » DEGRE T S • / / 3 <>  , 1H-,4X,  14HDR-LEOF  N FIJNCT  . 
111X,L7HPHI  - PARTIAL  SUM , 10X , 1HN, 4X , L4HLEGFNDRE  FUNCT, 

1 l IX , 17HTHFT A-PART I AL  SUM) 

J.,1  FHR.VAT  ( 4HII.  v=  , II  , 5H  MCP=,I1,4H  M1=,I2,5H  MVM=,I2,5H  MVR-,12, 

14H  KA=,F6.2,8H  T HF T AC = , Ft . 2 , 6H  NPHI=,I3,5H  NTH=,I3,4H  LX=,14> 

,02  FUKMATKH  LM=,Il,5H  MCR=»  I 1 » 4H  M1=,12,5H  MVM=,I2,5H  MVR=,I2, 

14H  KA=,F6.?.RN  THE T A C = , F 6 . 2 , 6H  MPHI=,I3,5H  NTH=,I3,4H  LX=,I4) 

, ;14  FORMAT  ( 10F  7. 2 ) 

FHRMAT (2F 15.7 ,F9.3»2L15.7*F9.3) 

, rb  F 0 RM  A T (//11X«9HFCT  ( IT  , L ) » 1 4X  , 2HI)8,  1 S X , 9HECF  ( 1 T , L ) , 1 4X  , 2HDB  ) 

, q7  FTP'- AT  ( / /I  1 X ,9  HER  \ ( I T ,L  ) , l 4X  * ?HPR  , 15  X,9HFRF{  1 T , t ) , 14X ,2HDR  ) 

, fjfc  F 'Rv  AT  ( ?X  , • ***THI  S PATTERN  CUT  IS  ZERO  AND  IS  NUT  PLOTTED.****  ) 
,«q  FORMAT ( IH1 ,30X t 'RADI  AT  ION  PATTERN*  ' , F 4 . T , 'DEGREE  CONE'// 

12X,  I1HEXC  ITAT10N.25HFLEM.  TYPE  ( C I R C = 1 * R A D=  2 ) = » I3»r>X, 

1 L6HE-C DMP(TH  = 1 ,PH=2)=»  13/1 3( 4X ,4HK A = , F 7 . 2 , 2X , 2HM=  , 13))) 
ini  FORMAT ( ?X,2F 10. 3 ) 

,12  FORMAT ( 2X»  ?I 3, 5HTEST  = ,E  fl.2 , 3X, 5HTFSF  = ,F8. 2) 

PIO?  = PI/2. 

STOP!  = DS0RTI2./PI) 
us  topi  = . 5 * s T np i 
TPI  = ?.*PI 
DTR  = P I / 1 80. 

IT  - 0 
CALL  PLOT  I 

1 RF AD ( 5 , 1 l 0,  END  = 99 ) THET A0,KA , THE  TA I , THt  T AF ,OTH, PHI , DPMI , DK , 
1M,MCF,MVM,MVR 

IF (THETA U.FO.U.. AND. KA.fcO.O.)  STOP 
IF (THETAO.EO.D.)  THE TAO  = THE  T AF 
WRI Tf ( 6, 109 ) 

HR  I TF  ( A, 1 10)  THE TAO, K A . THE T A I , T HE  T A F , OT H,PH[ ,DPHI »DK, MC  R , MVM , WVR 
IF  (DPMI  .EO.O. ) DPMI  = 5 . 

Ml  - M 
T T = T T + 1 

K Z R = 0.0 

NPHI = P H I /DPH I ♦ 1 . 

KAS I I T ) - KA 

ST  p I K A = SORT  ( TP  I *K A ) 

ST  POT  P = STP1KA/TPI 

STROP  = STPIKA/Pl 

SNTt  = DSIM  DTR*THET  AO) 

S\ TO  1 = MS TCP  I /SNTP 
SNT02  = STOP  I /SN'TO  * 

AN  SO  = K A * S N T O 
T AKSn  = 2.*AKS0 
Pf) DK  = P I /DK 
IM  = 0 

300  GO  T0( 310, 310, 330) ,MCR 
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MM  ? S T HE  NUMBER  Of  TERMS  CONSIDERED  IN  THE  M-SERIES,  AND 
IS  OFT  INED  PY  THE  OATA  VALUES  OF  MVM.MVP. 

,10  MM  = M V m 
G'l  TO  3 
,30  MM  = »'VR 

3 RM  = m 

I M = I M + l 

MS  1 I M ) = M 

EC  = POP  K * T AK  S 0*C0S ( R M * DK / T AKSO)/ ( (PODK*AKSO) * *2  - RM*RM) 
EC  l = T C * ST  POT  P 
EC?  = E C * S T POP 
F ( I T , IM ) = EC 

4 MO  = 1 

THETA  = T I ' f T A I 
IA  = 1 

r CALCULATE  ANO  STORE  COS  ( MPH  I ( 1)  ) , S I N ( MPH  I (I)  ) 

DO  2 l = 1 , NPH I 
R I = I 

PH  IN ( I ) = (PI  - 1 . ) *DPH! 

RAOVPH  = P.M*DTR*PHIN{  II 
SUP  U)  = SI Nl kADMPH) 

2 C M P ( I ) = COS { R ADMPH) 

C MX  IS  N° . OF  TERMS  IN  NU-SERIES  REU.  FOR  CONVERGENCE. 

MX  = KA  ♦ 2.*PM  ♦ 10. 

MMTF(6,115)  T HE  T AO » KA , M 

MP 1 = M ♦ 1 

ONU  = l 80 . Z THE  T AO 

DO  11  MD  = 1,2 

IF ( MD.EO. 1 . AMD.MCR.E0.3)  GO  TO  II 

5 XNI)  = PM  + .5  - ONO 
ID  r MD  - 1 

c CALCULATION  of  ZEROES  OF  LEGENDRE  FIJNGT  I HNS  , NU  , NUPR  I MP  . 

DO  10  J = 1 ,MX 

N'll  = XNU  DNU  - .04 

N<J2  = Mil  <■  .08 

CALL  ROOTS! T 1 1 FT  AO  » NU 1 » NU2 , M , XNU , ID) 

C ID  IS  PARAMETER  F PR  MODE. 

CALI  RFS IDUITHET AO, XNU, M, 10,RES( J,MM  I 
C EVALUATE  NU-COEf  IN  NU-SFRIES,  NOTING  THAT  RESIDI'E  FOR 

C RADIAL  CASF  IS  SAME  AS  FOR  MODF  2. 

TTD  = XNtJ/RE  S ( J,  MD  ) 

T = ITD/IXN U**2  - .25) 

TT  = T T n 

XN 0 I m = XNU*  P I 02 
X N P = r mrl X ( 0.  , XNP I M ) 

XEN  = CFXP(XNP) 

PR  (J,MD)  = T * XE  N 
IF (MD.NE .2)  GO  TO  10 
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(J,MD*l)  = TT#XEN 

,j  NtMJ.MD)  = XNM 
,"l  CONTINUE 

WPITF(6,118)  ( ( J, “D, NIK  J ,MD) , MO=l , 2 ) , J= 1 , MX  ) 

CALCULATION  OF  BESSEL  DEPENDENCE  TIMES  F { NU  ) /R  ES  I DUE  . 
NRITF(6»112)  mcr 
(in  16  NO  = 1,3 

(f  (ND.LF.2.AND.MCP.EQ.3)  GO  TO  16 
(F {MD.E0.3.AND.MCR.E0.2)  GO  TO  16 
WRITE(6,116)  MD 
K = M0 

|F (MU.EQ.3)  K = 2 

EXECUTION  OF  ST  500  INITIALIZES  AR(J,*D)  TO  PREVENT  OVERFLOW. 
IF (MD.FO.2.AND.M.EO.0)  GO  TO  500 
no  15  J = 1 ,MX 
XNU  = NUlJ.K) 

IF  (MLr.NE.3)  GO  TO  13 

CALL  BTSSIlKA, DK , XNU , B J ( J ) ) 

BISSI  EVAIUATFS  INTEGRAL  FOR  RAOIAL  SLOT  CASE 
GO  TO  14 

CALI  BESSEL  IN A, XNU , R J ( J ) ) 

IFI^D.FQ.ll  GO  TO  14 
XN'IPl  = XNU  *■  1. 

CALI  PESSEL ( KA,XNUP1 , B1  ) 

FOLLOWING  EXPRESSION  IS  FUR  D/DX  I X * JNl)  I X ) ) . 

RJ(J)  = (XNU  + .5J/KA  *B J ( J ) - B1 
1<,  IF  (J.LT.6.0R.J.GT  .MX -6 ) 

l WR  ITE  (6 , 130)  J,NII(J,K),RR(J,MD),RES(J,K),RJ|J) 

HJSIMG  = HJ ( J ) 

AR I J , MD ) = BR I J, MD) *RJS ING 

15  CONTINUE 
GO  TO  16 

600  CONTINUE 

DO  17  J = 1 , * X 

17  A«(J,MD)  = ZERO 

16  CONTINUE 

CALCULATION  F LEGENDRE  FUNCTIONS  AND  PARTIAL  SUMS. 

I A T = l.*  (THFTAF  - THETAll/DTH 

NAI  = TUETA1/DTH 
If  (NAI  .NE . I ) I A T = I AT  ♦ 1 
I AM  = 1 AT/2 

LOOP  1 r.  INITIAIIZES  the  nu-sepies  total  sum  fop  fach  new  m. 

nn  lb  I 3 1 , I at 

ACT(I)  = ZERO 

ACF(I)  = Z E RO 

AR  T ( I ) = ZERO 

18  ARE ( I ) = Z F RO 
CF  = RM/SNTO 
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DO  POO  MD  = 1,2 

If  ( wo.ra.  1 . ANO.MCR.EQ.3  » GO  TO  200 

in  wr  i if.  (*•  1 1 n uo, mcr 

TOFT  A = THfTAI 
DO  210  1 A = 1 , 1 AT 

GO  TO  (240,241),  MD 

’/,0  IF  ( 1A.FQ.1.0R.  IA.E0.  I AM  . OK  . I A . EO  . I AT  ) WRITF(6,120)  TMF  T A 

GO  TO  20 

,>41  1 F « I A.FO.  1 .OP  . 1 A.  EO.  I AM  .OR  . I A.  FQ.  I AT  ) WRITF(6,180)  THETA 

20  RAO  = DTR*TH1 TA 
SNT  = DS INC  RAD) 

CST  = DC  OS (RAD) 

ST ( I A)  = THE  T A 

THE  ST  1 , ST  2 , ST  3, ST  4 ARRAYS  AWF  CREATED  FOR  PLOTTING  ROUTINES. 

ST  1 ( I A ) = THETA 

ST  2 ( I A)  = THETA 
S T 3 ( IA)  = THfTA 
ST  4 ( I A I = THETA 
PS  = ZERO 
PSD  = ZERO 
RSDR  = ZERO 
PS  R = ZERT 
PSING  = 0. 

PDSING  = 0. 

p = o.no 

PI  = 0. DO 
PD  = 0 . DO 
DO  30  J = 1 , M X 
XNO  = NU(J.MD) 

CALL  l EGF. ND(  THFTA,XNU,MP1  , PD, 0,0) 

PDSING  = PD 
IF (”125,26,22 

22  CALL  LFGFN’DI  TOFT  A,XNO,M  ,P1 ,0,0  ) 

PD  =R”*CST/SNT *P1  f PD 
PDSING  = PD 

CALL  I EG END!  T HE T A , XNU , M , P , 0 , 0 ) 

PSING  = P 

If  (MCR.EQ.3)  GO  TO  26 
PS  - PS  ♦ AR ( J ,MD) *PSING 

25  PSD  = PSD  F AP ( J , MO  I *P0SING 

26  IF (md. NE. 2. OR. MCR. EQ. 21  GO  TO  27 
PSOR  = PSOR  ♦ AR ( J ,MD+ 1 ) *POS I NG 
PS  •*  - PSR  f AR  ( J.MDFl  )*PS  INC, 

PSOR  AND  PSR  ARE  TOTAL  SUMS  FOR  THF  RADIAL  CASE. 

27  IF  ( I A. TO. 1 .OR . I A . FO.  I AM.Op. I A. FQ. I AT  1 GO  TO  31 

GO  TO  30 

1 1 IF  I J.L T . 5. OR. J.GT.MX-8 ) GO  TO  32 
GO  TO  30 
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. > i f ( Mf.u  . fo.  3 ) on  rn  ?a 

rtKiTr(f>*i2i)  j.pd.pso, j,p,ps 

. |f  {MD.NI  .2. OR. MCR. FO. 2)  GO  TO  30 
WRITI  16,121)  j.po.psdr, j,p,psr 

,u  C'TNT  I NIIF 

F-  V ALU  AT  1 ON  OF  FIELD  COMPONENTS  AT  f AC  H THFTA  FOR  EACH 
C F 1 A R = WM/SNT 
C F 1 = CM  PI  XIO.  ,CF1  AR) 

GO  T0( 220,230) ,MD 


ACT!  IA) 

= 

ACT!  IA)F 

PSD 

ACE  ( I A) 

- 

ACF  ( IA>- 

PS*  FI 

GO  TO  215 

TO 

GO  TO (235 

,237,236) , 

MCR 

AC  T ( I A 1 

= 

AC  T ( I A ) «■ 

PS*CF 1 *CF 

acf  ( i a ) 

= 

ACF ( I A ) + 

PSD*CF 

3 (• 

AR  T ( I A ) 

= 

ART ( IA)  ♦ 

PSP  +CF  1 

ARE  ( I A ) 

= 

AR  F I l A ) + 

PSOR  *0J 

GO  TO  215 

Z ' 7 

AC  T ( I A ) 

= 

A C T ( IA)  t 

PS*Cf 1 *C  F 

ACF ( I A ) 

= 

ACF ( IA  ) + 

psd*cf 

; 1 5 

IF ( I A .FO 

• 

1 1 THETA  = 

FLOAT  I N A I ) * OT  H 

210 

THE  T A = 

THFTA  DTH 

2 10 

CONTINOF 

EVALUATION  of  AMPLITUDES  and  FINAL  FIELD  COMPONENTS. 

NTH-  I A T 

MT  = NTH 

LX  = M PH  IX NTH 

K Z R IS  A COUNTER  FOR  INITIALIZATION  OF  THE  TOTAL  FIELD  ARRAYS. 
K Z R = KZP  ♦ IT 
IFIKZP.GT.IT)  GO  TO  45 
DP  40  L = I , LX 
FC  T ( I T ,1  ) = ZERO 

FCFIIT.U  = ZFRO 
ERF  I IT ,L  ) = ZERO 
40  FRTI IT ,1 ) = ZERO 
45  CONTINOF 

I F I MCR .F  0.3 ) GO  TO  5h 

THE  Dn  5 G LOOP  EVALUATES  FIELD  COMP.  OF  CIRCUM.  SLOT. 
tCTM,  FCFM  APE  THE  FIELDS  AT  A POINT  FOP  A PARTICULAR  M. 

ECT,  FCF  ARE  THF  TOTAL  FIELDS  AT  A POINT,  CIRCUMF.  CASE. 

no  so  l = l , nphi 

no  51  J = 1 , NTH 
L = NTH*I  - NTH  + J 
IF(M.MF.O)  GO  Tn  54 
E C T M = rri*CVP(N*ACT(J) 

I C F ( I T , L ) = ( l .E-I5, 0. ) 

SF(L)  = -300. 

SFPIJ)  = SFID 
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Gn  T n 55 

. /.  f r iw  - i i j * ac  i ( j ) 

EOF*  = FC2*SMP( I ) *ACF< J ) 

l r { CA  li  S ( FCFM  ) .LT.  1 . E - 1 5 ) E C E M =(1.E-15,0.) 

ECF(  IT, l ) = E C F ( I T * L ) * F C F *4 

SFIL)  = 10  . * ALOG 101PEAL  (EC  FM*CONJ  Gl  ECF M ) ) ) 

SF  PI  J ) = SF  ( L ) 

r.c,  Ir  (CAPS(ECTM).LT  . 1 . E - 1 5 ) EC  r M =(l.E"15,0.) 

EC  T ( I T ,L ) = FCT(IT,L)  ♦ EC  T M 
S ( L > = 10.  *AL0G10 (RE AL { ECTM*CONJG( ECT u ) }) 

SP ( J ) = S ( L ) 

SI  CONTINUE 

USE  "WRITE (6, 170)  P H I N { I ) , { S P ( J ) , J =•  1 . NT  H ) " AND  FOR  " SF  P ( ) HERE. 
SO  CONT I NUE 


STATFMFMTS  THROUGH  *58*  ARE  ONLY  FOR  WRITING  OUT  y ODE  PATTERNS. 
WRITE! 6 *143)  1ST (J)  , J = 1 , ► T ) 

00  96  I = 1 , NPM I 
N I N = NTH* ( I - 1 ) 

N I N 1 = N I N ♦ 1 

N I N *■'  T = NIN  + MT 

So  W R l T E ( 6 , 170)  PHIN I I ) , ( S ( J ) . J = N INI ,N . NMT  ) 

IF ( M.EQ.O)  GO  TO  59 

WR I T F (6, 143 ) ( ST( J) , J=1 , MT) 

DO  5 8 1 = 1 * NPH I 
NIN  = NTH*  * 1 - 1) 

N I N 1 = NIN  + 1 
NINMT  = NIN  ♦ MT 

5R  '.RITFI6, 170)  P H I N I T ) , ( S F I J ) , J = N I N 1 , N 1 NMT  ) 


59  IF (MfR.EO.2)  GO  TP  69 

PRIM,  FRFM  ARE  THE  FIELDS  AT  A POINT  FOR  A PARTICULAR  m. 
FRT,  ERF  A R F THE  TOTAL  FIELDS  AT  A POINT,  RADIAL  CASE. 

DO  60  I = 1 , NPH I 

DO  61  J = 1 , NT  H 

l = NTH* I - NTH  + J 

IF  I M . NE . 0 ) GO  TO  64 

ERFM  = f M P ( I )*ARF (J)*SNT01 

ERTIIT.L)  = ( l.E-15,0.) 

SIL ) = -300. 

SP ( J)  = SIL) 

GO  TO  66 

6 4 ERF^  = C*P(  I ) *ARF I J) *SNTP? 

F R T M = SMP(  I ) *ART ( J ) *SNT02 

I F ( CAPS ( FRTM ) ,L T . 1 .E-15  ) ERTM  =(  l.E-15.0. ) 

FR  T l I T ,1  ) = I P T ( l T ,L  ) ♦ ERTM 

SIL)  = 10.*AL0G10(RF AL  (ERTM*CONJGIERTM)  ) ) 

SP  ( J ) = SIL  ) 


* 


6 b IF(CARS(EPFM).LT.1.E-15)  T P F M =(1.E-15,0.) 

F P F ( I T » L ) - FPF(IT,l)  * E R F M 

Sf  (U  = 10.*ALOG10(REAL (ERFM*CONJGIERFM) ) ) 

S F P ( J ) = SF(L) 

61  COiTIN'UF 

USE  "WPITF (6*170)  PHINl 1 ), (SP( J) ,J=1,NTH)  " AND  FOR  "SFP()  HERE. 

AO  CONTINUF 


STATEMENTS  THROUGH  *68*  ARE  ONLY  FOR  WRITING  OUT  MODE  PATTERNS. 
IF ( M.FQ.O)  GO  TO  67 
WR I T F ( 6 * 143) (ST(J)*J=1,MT) 

Of)  66  I = l.NPHI 
N I N = NTH* ( I - 1 ) 

NIN1  = NIN  4 1 
N I N M T = NIN  * MT 

66  WRITE  (6*170)  PH  I N I I ) , I S ( J ) , J = N I N 1 , N I NMT ) 

67  WRITE(6,143>  I ST ( J ) , J = 1 , MT ) 

no  68  I = 1 , NPH I 

NIN  = NTH* ( I - 1 ) 

N I N 1 = NIN  * 1 
N I NMT  = NIN  * MT 

68  WRITFI6.170)  PH I N ( I ) * ( S F ( J ) , J =N I N 1 , N I NMT ) 


69  IF(M.GF.MM)  GO  TO  ISO 
M = M 4 1 

GO  TO  3 

3b0  IF l VC  P . NE  . 1 ) GO  TO  360 
MCR  = 3 
ISW  = 1 

IF ( M.GE. ywR ) GO  TO  360 

M = M 4 1 

GO  TO  300 
360  CONTINUE 

C THE  "IT"  CASE  TOTAL  PATTFRNS  ARE  COMPUTED  AND  PLOTTED  HFLOW. 

WRITE (6, 112)  MCR 
T H T T A C = 180.  - THETAO 
IF  ( ISW.FQ.  1 ) Mf.R  = 1 

C THE  WRITE(IO)  STATEMENTS  WRITE  THE  COMPUTED  DATA  ON  TAPE. 

WRITE!  10)  MCR 
DO  A 00  LM  = 1,2 
M7  = m - m 1 + 1 

IF  ( I.M.  FC.2)  M z = MVR  - Ml  + 1 
WM  ff  (ft,  11?)  MCR 

If  (MfR.EO.?. AND.LM.EQ.2)  GO  TO  400 
IF  ( 4'C.R  .EC.  3.  AND. I M.EQ.  1 ) GO  TO  400 
KCR  = 0 

TO  370  I = 1 * NPH I 
Tf  ST  = -?H0. 
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TEST  = -200. 

nn  375  j = 1 , nth 

L = N T H ♦ I - NTH  + J 
GO  T0( 371 , 37? ) ,LM 

371  S { L ) =10.*ALDG10(PEAL (FCT  ( IT  ,L  )*CONJG(FCT ( IT  ,L  )>)  ) 

S F ( L ) =10.*AIOG10(RFAL(ECF(  I T » L ) *C ON JG ( EC  F ( I T « L ) ) ) 1 

GO  TO  3 74 

3 72  S ( l ) = 10. *ALOG  LOIRE  AL  ( E R T ( IT  , L ) *CON J G ( FP  T ( IT  fL  ) ) ) ) 

SF  (1  ) = 10.*A(_0G10(RFAL(FRF(  IT  »L  )*CONJC>(ERF(  IT*!.  ) ) ) ) 
3/A  IF ( S ( L ) .GT. TEST)  TFST=S(L) 

IF  ( SF ( L ) .GT. TE SF ) T E S F = S F (L) 

3 76  CONTINUE 

SM  T ( I ) = TFST 
S'-'FII)  = TESF 

W° I T F ( fc  * 1 92  ) LM, I , TEST , TESF 
370  CONT INUE 

GO  TO  (37 6, 377), LM 

376  WRITE (6*181 ) L M, MC R , Ml , MVM f MVR , K A , TH E T AC , NP H I f NT H, L X 
WP I TF  I 10  ) LM , NCR , Ml , MVM , MVR , KA , THE T AC , NPHI , NTH, LX 
WR ! T E ( 1 0 ) ( PHI  Nil) , I =1 , NPHI ) 

WR ' TF ( 10)  1ST  (J)  *J=1 * N T H ) 

WR  i TF  ( 6,  1 86  ) 

WP ! TF (6, 185)  (EC.TI  IT,L),S(L)»ECFl  IT»L)*SF(L)tL-l»LX) 

WR I TE ( 10)  (ECT( I T , L ) »S(L) »ECFI IT  «L )«SF(L) » L = 1 * l X ) 

GO  TO  378 

3 77  WPITE(6,181)  L M , MC P , M l , M VM  , MVR , KA , THE T AC , NPH I , NTH , L X 
WR I TE 11 0 ) LM,MCR,M1,MVM,MVR,KA,THETAC,NPHI,NTH,LX 
WPITF(IO)  ( PHINl  I ) , 1 = 1, NPH I ) 

WRITF(IO)  ( ST ( J)  , J = 1 ,NTH  ) 

WR I TF ( 6, 187) 

WR I Tf 16 , 185)  (FRT(IT,L),S(L),ERFIIT,L),SF(L),L=1,LX) 
WRITE  110)  {FRT  ( IT  , L ) , S ( L ) »ERF ( IT  * L )»SF(L  ) *l  = l»LX  ) 

378  CONT I NUF 

C LOOP  A 1 0 FINOS  IARG'  ST  VALUE  IN  S(L>  AND  SF(L)  ARRAYS. 

SMT  N = -200. 

SNFF  = -200. 

no  A l 0 I = l, NPH l 

IF ( SVT ( 1 ) .GT.SMTN)  SMTN=SMT(I) 

410  IF ( SMF ( I ) .GT.SMFN)  SMFN=SMF(I) 

SMTT  = SMTN 
S^FE  = SMFN 

WR  I TF  (6, 1*>1  J SMTT.SMFF 
SNOP  m = SMT  T 

I F ( S«F F .GT .SNORM)  SNURM=SMFF 
WRIT!  (6,101)  SVTT, SNORM 
3 79  KC R = KCR  ♦ 1 

WRITE  (6, 113)  M C R * I M , KCR 
GO  TO  (381, 382), LM 
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1 VJ  ° ! T L ( ( ,U0)  THETAC*LM,KCP»  IKASI  IT  ) , M S C I ) » F ( IT  , I ) , I = 1 * MZ  ) 

r,tl  TO  3 B 3 

W»IT(-(6«  189)  THF TAC ,LM,KCR, I KAS( IT) , MSI l ) , 1=1 ,MZ ) 
j WRITE  (6,143)(STIJ),J=1,NTH) 

CTF  = -70. 
on  3r«0  1 = 1 , NPH 1 
N l N = NTH* ( I - 1 ) 

NI  N 1 = M N + l 
MNMT  = MN  + NTH 
If  (KCR.FG.2)  GC  TO  386 
on  380  J = I, NTH 
L = NT H*  I - NTH  *■  J 

SVT  m ,SMT( n AR  F PHI-CUT  NORMALIZATIONS.  SMTT,SMFF  ARC  FOR 
COMPONENT  Nf'RMAl  ZAT  I ON.  S NO  RM  = M AX  ( S M T T , S MF  F ) OVERALL  NORM  AL  I Z AT  . 
S ( L ) = S ( L ) - SNORM 

SF ( L ) =SF ( L ) - SNORM 

SP 1 J ) = SI  L ) 

SOP  I J)  = SF ( L > 

^0  CONTINUE 

*85  WR  1 TF  I 6»  1 70  ) PH  I N I I ) , I S I J ) , J=  N I N 1 , N I NM  T) 

THF  F nL  L OW I NG  CAROS  CONTROL  THE  NUMBER  OF  PLOTS  DESIRED. 
IFII.GT.l)  GO  TO  390 

IF ( SP ( ? ) .L T .CTF . AND.SP ( 3 ) .LT .CTF . AND.SPI6 ) .LT.CTF ) GO  TO  387 
CALL  SYMBOL! 1 . * 7 . 5 i . 14,36*0. »-l  ) 

CAl  L L FGNDI (L.3,7.5,.l4,0.,l,,=,,l,PHIN(I),0) 

CALL  C PLOT  3 ( 1 , ST ,9. ,0. ,20. , 14HTHET A , DEGREFS,l4,7.,-2a.,4„, 

1 2F,HPr’WF  P IN  DB,  THETA  C OMPONF  NT  , 28  , NTH  , SP  ) 

3 P 7 IF ( SEP (2 ) .LT  .CTF  . AND.SFP I 3) .LT .CTF . AND.SFPI6 ) .LT.CTF ) GO  TO  390 
CALL  SYMBOL  ( 1. ,7.5,. 14,36,0. ,-l) 

CALI  LEGND1  ( l.3»7.5*.14*0.»l  t * - 1 » 1 » P H I N I I )»0) 

CALL  C PLOT  3 ( 1 , ST , 9. ,0. , 20. , 14HTHFT A , DEGREES, 14,  7.»-28.  » 4 . * 
126HPOWER  IN  DB,  PHI  COMPONENT , 26 , NTH, S FP ) 

GO  TO  390 

i B 6 WRITFI6,170)  PHI N (I) , I S F ( J ) , J = N I N 1 , N I NM T ) 

boo  CONT  I Nt if 

IE  (KCP..EO.  I ) GO  TO  379 

THE  FOLLOWING  CARDS  PLOT  PATTERNS  Of  THETA  VS.  PHI. 


IF  I JS.EO.O)  GO  TO  40  0 
IF  ( JL .GT .NTH)  JL  = NTH 
0°  690  J = JS, JL » J ST  F P 
DO  6 PC  | = 1 , NPH I 

l = NTH* I - NTH  ♦ J 
SP ( I ) - S ( L ) 

SF  PI  I ) = SF  ( L ) 

880  CONTINUE 

IF (SPI2) .LT.CTF. AND. SPI 3). LT.CTF. AND. SPI6I.LT. CTF)  GO  TO  587 
CALI  SYMBOL (1., 7. 5, .14,37,0.,-!) 


CALI  Lf  C\('  1(1.3,  7.5,.14,0.,l»'  = 'tl,STIJI,0l 

CAlt  CPL  'T  3 ( l»PHIN,9.,0.,?0.*l  2HPH  I , D E GR  E F S * 12, 7 . , - 2 8 . , 4 . , 

1 2RMPOWE M IN  DP,  THETA  COMPONENT ,28,NPM1 ,SP) 

587  I»  ( SFP ( ? ) .L T . CTF . AND.SFP ( 3 ) .L T .CTF . AND. SEP (6  ) .IT .CTF ) GO  TD  590 
CALL  SYMPOL ( 1 . , 7.5, . 14, 37, 0. , - 1 ) 

CALL  LEGND1 I l.  3, 7. 5,.  1 4 ,0.  ,1,'  =•  , l.ST  (j|  ,01 

CALL  CPLGT3( 1 ,PHIN,9., 0. ,20. , 12HPH1 , DEG&FES* 12, 7. ,-28. ,4. , 

1 2 6HPPWE  R IN  DB,  PHI  COMPONENT , 26  « NPH l ,SFP| 

5 JO  CONTINUE 


400  CONTINUF 
GO  TP  1 
99  CnNT I NUE 
CNOFILE  10 

CALL  PLOT (0. ,0. ,999) 

STOP 

END 
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SUDr<  f • 1 1 I f-  I HI  G nI  [ ( X,  Nil,  J MIX) 
l)i'i)P|  [ PI-  I r IS  | IN  Nil , J Nil  X , XX 
PEAL*B  DBSJ1 
XX  = X 

IF  { X .r,T  . 100.  .E'R  .NU.GT  . 100.  ) r,(l  Tn  100 
JNIIX  - DRSJl  (XX, NO) 

RETURN 

100  CALL  PCSF  XT  (XX, NU,  JNlJX  ) 

PL  TUP N 
END 


SUBROUTINI  BE  SE  XT  I XX,  NU,  JMJXl  ) 
r.fiMMON/P.FSwKK/L  1 , L ?.  , A ( 700) 

REAL  * B XX,  Nil,  JNUXL  ,A 

THE  RFSFXT  SUBROUTINE  RETURNS  THE  VALUE  OF  THE  BESSEL 
c FUNCTION  OF  ORDER  P.LT.400  AND  ARGUMENT  X.LT.400.  LARGER  ORDERS 

C AND  ARGUMENTS  ARE  CONSIDERED  BY  EXTENDING  THE  DIMENSION  OF  A, LI. 

LI  = 700 
K - NO 

C A l L DB  S v/  2 ( X X , NU  , A ) 

JNUXL  = A ( K ♦ 1 ) 

Ri TURN 
EA'D 


S' MR  OUT  I N t Bf  RSI  (KA.DK.NUtMNT) 

IMIS  IS  AN  I N T F GW  A T I UN  PROGRAM  USING  SIMPSON'S  3 / R RULE. 
NR  IS  T OF  AIU*»fP  OF  INCREMENTS  ON  RK  AXIS,  MULTIPLE  OF  3. 
Dt  MF  NS  ION  f ( ) 

DO UP | F P R r f I S I UN  NU, «INT,  PJ 
RE  Al  K A 
DATA  NR/IN/ 

R 1 K = KA  - DK  / ? . 

R ? K = KA  ♦ OK/?. 

RN  = NR 

H»  = ( R 2 K - PlKj/RN 

A = 0.0 
PK  = R 1 K 

D)  30  I = If  NR,  3 
D ’ ?0  J = 1 , '♦  , 1 
k K P = P K * * l — 1 .S) 

CAl  l BTSSri  ( RK  , NU  , B J ) 

GnK  = COS  IRK  - KA) 

F ( J ) = R K P ♦ P J * l.OR 

l)  ( J.F  (.'.*)  GO  TO  30 
?0  Kn  = P K ♦ HR 

30  A - A ♦ I ?./Q.  )»M«*(f  ( I ) ♦ 3 . *E  ( ? ) ♦ 3 . *F  ( 3 ) ♦ F ( '♦  ) ) 

BINT  = A 
RE  TURN 
END 
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SU'iP'MII  1 * I l f . I MM  1 H>-  T t.  , MU,  H , PNU  , PI  , ! P ) 
i)  >ijri  F PRtCISION  P 1 , P2  i CF  , DA  , A , PNU  , R VGM 

:)iH,r'i 1 prfcisi  >\  pi  ,csoa  , snoa.csa,  SNA , t , x ,f  ,rad,  p,  c , tx  f y , no 

i ) \ T A PI/ 3 • l A 1 N C ? *>  3 5 8«) 7 9 / 

7 AP  = THfTA*PI/lRO. 

Y - ?.  *PS  I V ( P A l>) 

,)/>  : P Ml  ♦ 

J M r fj 

\ = (70.+PM-fP‘-’)  -\n 

IH  (NU.GT  . (69.  ♦•<«♦•<“)  ) A = -1 
SM  P A = l S I N ( 0 A | 
fS'iA  = PC  I S ( I A ) 

R\  = N 

Nil  = NH  ♦ (CN+2.  > 

J = 0 

r m 2 = k w * t-  M 

J f I Nl'OP  AO.  I T ) GO  TO  SO 
10  N'l  = NO  - I. 

\ = N(  K1  K A 1 1 ♦ ( .25*  . S*P*  ) *P  I 

SNA  = LS  1M  A ) 
fSA  = OCPS(A) 

C = 1. 
p?  = PI 
T = l. 

D1  r SNA 

70  T = -TM  (C-.S  ) *+2  - FV?)/(C*Y*(N'J  + C)  ) 

SNA  = SAAtG.SOA  ♦ CSA+SNOA 
rSA  = (CSA  - SNA*SNPA) /CSOA 
PI  PI  ♦ T * S N A 
C - C ♦ 1 . 

IF  ( PA  1S(T)  .GT  . l.n-o.  ANI1.C  .LT  . 35.  ) GO  TCI  20 
1 = . 7 r>  - f M 2 
n |\.FQ.(-1I)  GO  TP  AO 
J = J + 1 

IF  I J.tO. 1 ) GO  TO  10 
r.  - Nil  ♦ ?. 

TX  = 2.*OCOSlK  Af)> 

NT  = ( N ♦ 1 ) /? 
on  30  J = 1 , NT 
C = f - 2. 

P7  = TXtPl  - I 1 . 4-R/ (C*  ( C*  1 .)))  *P2 
TO  PI  = T X * P?  - ( I.*-R/(C*(C-l.)  ) )*P1 
J ? - M0l)  ( N ♦ 1 » 2 ) 

NU  - A'U  - «N 
IF  (J.tO. 2)  P 1 - P2 

IF(IP.FP.l)  Pf TOON  TO  SHPROOTS  FOR  COMPLETION  OF  ROOTS. 

AC  I p ( I P . F („  . I ) RF  TORN 
X - Nil  *RM  - . c> 
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T - Ml  +rm  _ #t, 

I M M . E Q . 0 ) X = 1 . 

IF  ( M.L  T .2)  GO  TO  4 5 
DO  44  I = 2,V 
T = T-l. 

4 4 X = X * T 

•♦5  IF  ( NO  . G T .110.)  GO  TO  46 
CALL  GAVPAT (NU.M.GR, I ) 

CF  = 2./DSQRT (PI * Y ) * X*GR 

PNU  = CF*P1 

RETURN 

46  T = l./( 16. *NU ) 

CF  = ?./nSQRT(PI*Y*NU)*X*(  l.-?.*T*(  l.-T)  I 
PNU  = C F * P 1 
PF  TURN 

60  NO  = NO  -RN  - ?. 

INU*PAD.LF .9 » NOW  EXECUTING. 

F = MJ*N'U  - .?  5 
X = DC OS  I F AD ) 

Y = ( 1.  - X) /2. 

PI  = 1. 

T = 1. 

C = 0. 

60  C = C *■  1. 

T = -TMF  - C*  (C-l  . ) ) / ( RM  + C ) *Y/C 
PI  = PI  + T 

If (OAPSIT). GT. 1.0-9. AND. C.LT. 25.)  GO  TO  60 
call  GAMRAT (NU,M,GR,2) 
rmgn  = rm  ♦ i. 

GM  1 = GANMAR  ( Rf>GM) 

IF (NU.LT. 100. ) CF  = GR/GM1 
IF (NU.GE.100. ) CF  = NU**(M+M) /GMl 

PNU  = C F * ( ( ( 1 . -X 1 / ( 1 ,+X I ) **( ,5*RM ) )* ( ( -1  . ) ) *P 1 

RF  TURN 
END 
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SimPO'JT  1NE  KFSIOUIT  'ETA,  NU  , M , lO.DPDNUl 

SUR-PES I DUE  RETURNS  DPDNU , THE  T AD . WHEN  ID=0  AND  ITS 

THETA  DERIVATIVE  WHEN  !D=1. 

OOURLE  PRECIS  ION  PI  , F ,R  AD  ,C  S T , SNT  * OPDNU  , DPNU20  ,Ntl , P2  , PN'J2 
DI  ‘TENSION  I (4  ) 

DATA  PI ,H/ 3. 1 59265358979, . 0001/ 

NU  = NU  - 2 . *H 
I = 1 

10  call  LEGEND! THETA, NU ,«,F (I) , P2,0) 

IF l I P.FO.O)  GO  TO  15 
M P 1 = P+1 

CALL  LEGEND! THETA, NU»WP1,PNU2,P2,0) 

RAD  = THET  A*  P I / 1 80. 

CST  = DCOS(RAD) 

SNT  = PSIN(RAD) 

RM  = M 

F ( I ) = PN*  CST/SNT*F(I)  + PNU2 
151=1+1 

NU  = NU  ♦ H 

IF  ( I .FO.  3)  NU  = NIJ  + H 
IF!  I ,LT  .5)  GO  TO  1C 

DPDMJ  = ! F ( 1 ) - 8 . «=  F ! 2 )+8.*F  ( 3 )-F  ( 4)  )/ ( 12. *H) 

NIJ  = N'l  - 3.*H 

RE  T UP  N 

END 


FUNCTION  GAMMAR(X) 

POURLE  PRECISION  X 
1FIX.GT.1.)  GO  TO  2 
Z = X 

1 IF ( l .GT . 0. ) GO  TP  3 
Z = Z + 1. 

GO  TP  1 

2 Z = X - 1. 

3 T 1 = ? +.5 
TZG  = T1  ♦ 5. 
n = T ZG  * * T 1 

T 1 = FXP ( -TZG)  *T l *2. 50662027465 

GAvz  = T l*( 1 . + 76. 18009173/lZ+l  . ) -86 . 50 5320 3 3/ I Z + 2 . ) 
1+24. 0140r)£,22/(  Z + 3.  ) -1.23173 95 16/IZ+4.  ) +.  I 2085800  3E-2/ 
l!Z+5.i  -.536382E-5/! Z + 6.  ) ) 

IF  (X.GT.  1.  ) GL)  TP  5 

4 GAMZ  = GA^Z/7 

IF ( Z .CP. X ) GO  TO  5 
1=1-  1. 

GO  TP  4 

•-  GAMVAP  = GAMZ 
RETURN 
ENP 
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S UUP  OUT  INF  RO0TSm.ETA,NUl,NU2,M,N'l,lD) 

ROnT  S nr  PNU  ARF  OHT  A I NED  WHf  N ID  = 0,  OF  DPDTHETA  WHFN  10=1 . 
OOiJPLt  PR  E C I S I ON  PI, EPS 

DOUR  IF  PRECISION  NU,NU1  ,NU2,NUL  , P 2 , 0 , PNU  , R AO,  EP  , FPS  , FPL  ,SCT  ,SNT 
DATA  P I , EPS/ 3. I 41  592  65  3 5 t 979,  1 .D-7/ 

26  FORMAT ( I 4H  I T N , NU , F P , E P S , I 5 , 3E  1 5 . 8 ) 

NU  = Nl  1 1 
ITN  = 0 

10  CALL  LEGEND! THETA, NU,M, PNU, FP,  1) 

If  ( I D. E0.0)  GO  TO  15 
M P I = M 1 

CALL  L F GE  NO ( THETA, NU , MP I , PNU.P2, l ) 

RAD  - THET A*PI /180. 

CM  = DCOSIRAO) 

SNT  = DSIN(RAD) 

RM  = M 

FP  = RM^CST^FP  ♦ (NU+RM+-.6)  *SN  T * P 2 
15  0 = MJ2 

ITN  - ITN  * l 

IF ( ITN. EQ. I ) GO  TO  20 

FPS  = FP  - FPL 

IF IDA3S! FPS) .LT.EPS)  FPS  = DS I GN ( E PS , F P S ) 

Q = NUL  - (NU  - NUL)/FPS*FPL 
20  NUL  = NU 
FPL  = TP 
NU  = 0 

IF  (OARS!  NIJ-NUL). LT.EPS)  RETURN 
IF ( ITN.GT.20)  GO  TO  30 

GO  TO  10 

30  WR ITF (6,26)  I TN , NU , r P , F PS 
RETURN 
END 
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S U M P U U I I N f GAMRAT I NU, M, GR , IG) 

SUO-GAMRAT  CALCULATES  AND  RETURNS  THE  RATIO  OP  TWO  GAMMA 
FUNCTIONS..  C,  ( NU  + • 5 ) / G ( N U + 1 . ) AND  G ( NU  *RM  4 . 5 ) /G  ( NU- R M 4 . 5 ) . 
DOUHLF  PRECISION  NU,X,D1 
10  FORMAT ( 3E?0. 10) 

RM  s M 

IF  I IG.FQ.2)  GO  TO  30 
KNU  = NU 

D = NU  - FLOAT (KNU) 

01  = D + 1. 

X = 0 4 .5 
RATI  = X/Dl 

Y = GAMMAR ! X) /GAMMAR  C D 1 ) 

If  l KNU.  F 0.0)  R AT  l = 1. 

Y ^ Y*P  A T 1 

IF (KNU.LT. 2)  GO  TO  25 
DO  20  I = 2, KNU 

m = n i 4 1. 

x = X 4 i . 

PAT  = X/Dl 
20  Y = Y*RAT 
2 5 GR  = Y 
RETURN 

30  IF(m.EO.O)  GO  TO  50 
S = NIJ  4 RM  - .5 

T - S 

MT  = 2*M  - 1 
00  AO  I = 1,MT 
T = T - 1. 

AO  S = S*T 
GO  TO  60 
50  S = I. 

60  GR  = S 
RFTUPN 
END 
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appendix  b 

COMPUTER  PROGRAM  FOR  MODAL 
OR  TOTAL  FIELD  PATTERNS 

This  program  is  also  written  for  the  IBM  370  system  and  provides 
complete  field  patterns  and  plots  selective  cuts  of  modal  or  sum  patterns. 

The  field  component  data  is  provided  by  one  to  three  magnetic  tapes 
generated  by  the  general  program  of  Appendix  A.  This  auxiliary  program 
allows  for  detailed  study  of  the  modal  series  without  costly  recomputa- 
tions  of  the  general  program. 

Most  of  the  program  parameters  have  the  same  meaning  as  in  the 
program  of  Appendix  A.  Other  parameter  definitions  are  as  follows: 

NTAPE  = number  of  data  tapes  to  be  read. 

If  NTAPE  = 1 only  tape  identified  by  NC 
will  be  read. 

NC  = Tape  unit  number. 

IT  = Counter  for  the  total  number  of  complete 

files  read  from  all  the  tapes.  Must  be 
< to  first  subscript  of  the  E-arrays. 

Scaling  and  labeling  information  must  be  provided  in  the  arguments  of  the 
plot  subroutines.  Multiple  graphs  will  be  drawn  with  the  present  listing  in 
all  plots.  The  complete  listing  of  the  program  is  included  in  the  following 
pages. 


, I'tWfR  PATTERNS  AND  PUTS  OF  CONICAL  AkK/.Y  ^ k Li  l-  L t M « 

C JlMtur  If.  T T . 1 PROGRAM  IS  TAKEN  FROM  TARE  . 

COUPLE  X t 07  t 1 i , (.  6 5 ) » E f.F  ( 1 1 ? 665  ) , ERT  (11  ,6r5  ) , t RF ( 1 1 , oo5 ) 

L 0 M P L X foil  , ELEF,fckT(  , F P F F 

f<  t A L K A • iC  A A 

( OMMGN  /Lt-Lu/ECTT  Uo5)  , Ff  FF  ( 06  5 ) /RFLO/E  K T T ( OCj  ) , F R F F ( c 6 6 ) 

C.OfIMuN  /PVS/  PH1N(2G  ) ,ST(40)  ,S(665  ) ,SF  (t  l 3) 

COMM(iN/KANA/KAA(  13) ,M1A (15) ,MVMA(  15) ,KVkA  ( 15  ) 

PATA  NTAPf/i./ 

0 A i A i^TAPt/1/ 

1«1  PORMATl  AMI  LM=  , 1 1 ,5H  Mf.R=,!l»AH  M1=,I?,dH  MVM-,1,,5H  MVR=»I2, 

IaH  KA-.F6.2,  c.H  IHLTAC=,F6.2,6H  NPHI=,I3,5M  MH=,I3,4H  LX=,14> 

1 ■(/  FCRHATIAH  LM-,  II  , 3 H MCK=,T1,AW  Ml  = , 12, 3H  MVM=,12»5H  MVF  = * I 2 * 

1AM  KA=,F£.2,8H  f HE.  T AC  = » F 6 • 2 » 6H  M'MI - , I 3 , 30  NT)!=  , 13  *««H  LX=,IM 
1 6z  Format  ( l i o . /,fv,.3»2E13.7,P'7.3) 

l?c  FORMa  I ( n\  1 X ,YnFcT(  1T*L  ) , 14X,  7 HPh  ,1vX,F|'i  LHIT.U,  14X,2HD6  > 

1<<  7 FORMAT  l //l  IX  ,9Ii-kT(  IT  ,L)  , l*Xf  ?H05»  1 3X  , 9P,.  kE  ( IT ,L) , 14X,2h0F.) 

1 SR  FORMA  I t‘*F  1 r . 7 ) 

CALL  PL0T1 
k t - If 
IT  ^ 0. 

1 A 1 - I ( ♦ i 

1 F < IT  . f C . 1 1 . ANt  . M . F Q . 1 1 ) 00  10  3 1 C 
RL  AG  ( NO  t FUG-3  10  ) NOR 

Of  TO  ( j.j  1 ,301,302)  ,MCR 

C il  i~AN  Ri  F k '2  U.  NT  A MODE  OR  A DIFFERENT  LL'MtNT  POSITION. 

301  ( AO  (NO  ) lM  F. , Ml  , MVM  , MVR  ,KA,THETAC,  NPF  I , NTH  , LX 

■<i  AO  (Tic  ) ( FniM  1 ) , 1 = 1 ,NPHI  ) 

At  AO (NO)  < 51  ( J ) ,o-  1 ,NTH) 

f:t  AiMiiC)  (i.  L((II,L),S(L)fFCF(lT,L)»r>F(L),L-l»LX) 

P A A ( IT ) - KA 

MAUI)  - M 
MVMA(IT)  - K V M 
mvp.a(ii)  = r:vr< 

. U5‘  -VOll  (<.,le>)  ( lC  I ( 1T,L)  ,S(L)  ,t  ' E ( IT  ,L  ) ,Sr(  L)  ,L  = 1 ,LX  )"  HLRt  . 

I M.'lCR.  EG.  1 > GO  TO  30,. 

&n  T'l  3 CO 

3 02  READING  ) LMf  MCk  ,M1  , MVM,  MVR.KA  , TITET  Au,NPMI  , NTH,LX 

Rt  ad  om ) ( y y in(  i ) ,i  - 1 ,nf)ji  ) 

Rt AL (N(  > ( ST ( J ) , i ,NTH ) 

AD  ( NC  ) |iK1(II,L),S(L),0»F(IT,L),SF(L),LM,lX) 

V.  A / ( 1 T ) - K A 

V.  1 A,  ( IT  ) --  Ml 

MVMA(.T)  - MV  A: 

MV°  A ( 1 T ) - MVk 

c IG,  -WRI1.  (r»,U'5)  ( LfcT ( IT ,L ) , S(L) ,rkF ( I T ,L ) ,Sr( L ) ,L-1 ,LX)"  HFRr  . 

300  CO  TC  1 

3 1 u 11  - il  - 1 
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IF ( NTAPL.EQ • 1 ) GO  TO  312 
IF(NC.EQ.l?)  GO  TO  31? 

IF(NC.EC.ll)  GO  TO  311 
MC  = 11 
GO  TO  l 
,11  KC  = 12 
GO  TO  1 
q?  CONTINUE 

CO  390  II  = I » I T 

I F ( MCR  * EG . 3 ) GO  TO  317 

DU  316  I = 11,11 

DO  315  L = 1 » LX 

ECTT(l)  = ECTT(L)  ♦ ECT(I,L) 

-15  FCFF(L)  = c C F F ( L ‘ ♦ ECF(I,L> 
jl6  CONTINUE 

317  I F ( MCR . EQ. 2 ) GO  TO  330 
DO  321  I = II, II 
DO  320  L = 1 , LX 
ERTT(L)  = ERTT(L)  ♦ ERT  f I » L ) 

120  ERFF(L)  = ERFF(L)  ♦ ERFCI.L) 

321  CONTINUE 
330  CONTINUE 

THE  FOLLOWING  "IF"  STATEMENT  IS  TO  BE  UStD  ONLY  WHIN 
HIGHER  ORDER  MODE  PATTERNS  IN  SUCCESSION  ARE  REQUIRE!  . 
CALL  PHICUT ( Il,MCR,Ml,MVM,MVR,KA,THET AC.NEHI , NTH, LX) 
IF(II.EQ.3.0R.II.EQ.6.0R.II.EQ.9)  GO  TO  390 
CALL  P L OT 1 —37 « , 0 . , —3 ) 

390  CONTINUE 

CALL  PLUT(0.,0.,999) 

999  CONTINUE 
STOP 
END 
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SUBROUTINE  PHlCUt ( I T , IXR « M l .MVM , MVR , K A , THE T AC ,NPHI , NTH, LX ) 
DIMENSION  SMT(20) ,SMFf 20) ,SP( AO) ,SFP(40) 

COMPLEX  tClT,ECFF,LRTT,ERFF 
REAL  KA.KAA 

COMMON  /CFLD/LCTT (665 ) , FCFF ( 665 ) /RFLD/E RTT ( 665 ) , FR  FF  ( fcfcS  ) 

COM.:  JN  /PVS/  PHI N (20)  , ST( 40 ) ,S(665)  , SF  (665  ) 

CC  ■ iON/KAM A/KAA (15) ,H 1 A ( 1 5 ) ,NVH A ( 15) » MVR  A ( 15 ) 

...EXECUTE  'DATA  JS, LL , JSTEP/O ,0 , 1 /•  LAST  IF  THETA-CUTS  NOT  DESIRt 
DATA  JS  * JL  » J STE  P/  0,  0,1/ 

OATA  JS,JL,JSTEP/  7,17,1/ 

1 j3  FORMAT (2X, 1 2HMCR , LMP , KC R= , 3 I 5 ) 

j-,0  FORMAT  ( 1H1  »30X»  'RADIATION  PATTERN,  ' , F A . G , 'DECREE  CONE'// 

12X.11H EXCITATION, 2 SHE  LEM.  TYPE ( C 1 RC  = 1 , RAD=2 ) - , I 3 , 5X , 

1 1 SHE -COMP ( TH  = 1 ,PH  = 2 ) = , I 3/ ( 3 ( 4 X , 4HK A =,F7.2,2X,3HMl=,l_,c/t 
14HMVM=,13,2X, AH MVR  = , I 3 ) ) ) 

1^3  FORMAT  (//IX,  1CHPHI,THFTA=,F5.1  ,16F7.1/(9X,17F7.1)) 

170  FORMAT! Ffa . 1 , 3X , 17 F7 . 1 / ( 9X , 1 7F7 . 1 ) ) 

]M  FORMAT ( 4H 1 L M = , I 1 ,5H  MCR=,I1,4H  M1=,I?,5H  MVM=,I2,5H  MVR= , i ' , 

1 AH  KA= , F6 . 2 ,6H  TH E T AC = , F 6 . 2 , 6H  NPHI=,I3,5H  NTH=,I3,*H  LX-,1-) 
1:0  FORMAT (2El5.7,FY.3,2E15.7,F9.3) 

191  FORMAT ( 2X*2F10.3 ) 

192  F0RMAT(2X,2I3,5HTEST=,F8.2,3X,5HTESF=,F8.2 ) 

T 2 FOLLOWING  TWO  STATtMENTS  DETERMINE  CURVE  IDENTIFIER 

.....FOR  THE  »LSEG»  PLUTT1NG  R 01  IT  I NF 

I T 1 = IT  - 1 
KL  = MOO ( I T 1 , 3 I 
DO  400  LMP  = 1,2 

IF (MCR. EL. 2. AND. LMP. EL. 2)  GO  TO  400 
IF (MCR.fcf.3. AND.' MP.EC. 1 ) GO  7 C 400 
KCR  = 0 

DU  370  1 = 1 ,NPH 1 
TEST  = -200. 

T F.  SF  - — <;uG  • 

DO  375  J = 1 , NTH 
L = NTH* I - NTH  ♦ J 
GO  T0( 371,37?) ,lMP 

?71  S ( L ) = 10.*ALUol0(REAL( ECTT(L) *CONJG (fcCTT(L)))) 

S F ( L ) = 10 . * AL0G1 0 ( R E AL ( ECF  F ( L ) *C0N JG ( ECFF  (L)))  ) 

GO  TO  374 

372  S(L)  = 10.*ALUG10(REAL(EPTT (L)*CONJG(ERTT (L ) ) ) ) 

SF ( L ) = 10.*ALUG10(REAL(ERFP(L )*C0NJG(ERFF(L  ) ))  ) 

374  1 F(S(L ) .GT .TEST ) TEST=S(l) 

IF(SF(L).GT.TESF)  TESF^SF(L) 

375  CONTINUE 

SMT ( 1 ) = TEST 
SMF(I)  = TE  SF 

WRITE(6,192)  LMP, 1, TEST, EESF 
370  CUNTINUE 
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WRITE<6, l<;i ) (SHT(l),SMP(n,l  = l,Nr  Mj  ; 

00  TO  ( 376 ,377) , IMP 

.1/0  WK  II  L ( 6 , lb  1 ) LHr  ,MCR»M1A(  IT  ) ,MVMA(I  ( ) , H V* , ( . , ) , n- , ( i i ) , 

1THET AC , NFH 1 ,K1H,LX 

USE  "WRITE  (6,183)  (ECTT(L)»S(L  ) ,ECFF(  L)  *SF  ( L)  , L - 1 » L X ) " ML*..  . 

GO  TO  376 

377  V/RITE(6,18  I ) LMP  , NCR  , Ml  A { IT  ) » HVMA  ( IT  ) »MVR  A ( IT  ) ,KAA  ( II  ) , 
1THETAC.NPHI ,NTH,LX 

USE  "WRITE(6 , 165  ) ( E RTT  ( L ) , S ( L ) , ER  FF  ( l ) , SF  ( L ) , L = l , L X ) ••  HER  t . 

378  CONTINUE 

LOOP  410  FINDS  LARGEST  VALUF  IN  S(!  ) AND  SF(L)  ARRAYS. 

SMTN  = -?0(  . 

SMFN  = -2 00. 

LU  410  I = I , Nz'hl 
IF(SMTd).uT.SMTN)  SKIN=SMT(I) 

410  IF(SMFU).C-T.SMFN)  SMFN=SMF(I) 

SMTT  = SMTN 
SMFF  = SMFN 

WF.  Hi  (6,141)  SMTT,  SMFF 
SNOKM  = SMTT 

IF( SMFF .GT . SNQKM)  SNORM=SMFF 
WRITE(6,191)  SMTT,SNORM 
370  KCR  = KCR  ♦ 1 

WRITfc(6,113)  MCR,LMP,KCR 

WRITE (fc, 140)  THLT  AC  »LMP  »KCR , ( K AAl I ) , M 1 A ( 1 ) »MVMA (1 ) ,MVRACI) ,1  = 1,  IT) 
WRITE (6, 14 3)  <5T< J ) ,J=1,NTH> 

CTF  = -70. 

DO  390  I = 1 , NPF I 
NIN  = NTH* ( I - 1 ) 

N IN  1 = NIN  ♦ 1 

NINMT  = NIN  ♦ NTH 

1 F < KCR . EG • ? ) GO  TO  386 

"SNURM"  NORMALIZES  PATTERNS  WITH  RESPECT  TO  LARGEST  VALUE 
f.  F BOTH  FIELD  COMPONENTS.  "SMTT"  AND  "SMFF"  ARE  FOR 
I NG 1 V I DEAL  COMPONENT  NORMALIZATION. 

CO  360  J = 1 , NTh 
L = NTH*  I - NTH  + J 
S I L ) = S ( L ) - SNURM 
SF  ( L ) =SF(L)  - SNORM 
SP( J)  = S C L ) 

SF  p ( J ) = SF  < L ) 

• C f NT INUE 

w 1)1(6,170)  PHiN( I ) , ( S( J), J=NIN1, NINMT) 

( I .GT . 1 . AND. i .LT.NPHI  ) GO  TO  390 
: ISPI2).LT.CTF.AND.SP(3).LT.CTc.AND.SP(6).L(.CrF)  GO  TO  387 

It  ,*  : ) .L  T .CTF.  AND.SFP<  3)  .LT  . CT  F . AND  . SFP  ( 6 ) .LT.CTF)  GO  TO  391 

TO  190 

*,,1/0)  PHIN(  I ) , (SF(  J)  , J*NIN1  , NINMT  ) 
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3c;0  CONTINUE 

IF(KCR.EQ.l)  GO  TO  il<) 

CALL  SUB-FBCUTIKL, ) FOR  CONTINUOUS  PLOT  OF  PHI=0  ANL  PHl=lbO. 

CALL  FDCUT(KL,NTH,NPHI,LX,LHPI 

THE  FOLLOWING  CARDS  PLOT  PATTERNS  OF  THETA  VS.  PHI. 


IF(JS.EQ.O)  GO  iJ  591 
I F ( JL.GT.NTH)  ’L«NTH 
DO  590  J » JS,  'L,JSTEP 
DO  580  I « ltNrHI 
L = NTH* I - NTH  ♦ J 
SP 1 1 1 = S(L) 

SFP(I)  • SF ( L ) 

«,fl0  CONTINUE 

CALL  PBPL0T<KL,PHIN,9.,0.,20. ,12HPHI , DEGREES, 12, 
?SP,7.,28.,-4.,28HP0WER  IN  DB,  THETA  COMPON  ENT  , 2 8 ,N  PH  I , 1) 
CALL  PBPLOT (KL,PHIN,9.,0.,20. ,12HPHI»  DEGREES, 12, 

1SFP,7.  ,28.  , —4  • .26HP0WER  IN  DB  , PHI  COMPONENT, 26, NPHI  ,1  ) 
IF ( IT .E0.3.0R. IT .E0.6.0R.IT .E0.9)  GO  TO  590 
I F ( LMP . EQ . 2 ) CALL  PLOT ( -24 . , 0 . ,-3 I 
CONTINUE 
CONTINUE 


I F ( IT.EQ.3.0R.IT.EQ.6.0R.IT.EQ.9)  GO  TO  400 
I F ( LMP . EQ . 2 ) CALL  P LOT ( -1 3 . , 0 . , -3 ) 


40C  CONTINUE 
900  RETURN 
END 
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SUBROUTINE  PBPLOT ( K , X , XN , XS , DX , XL AB , NLX , Y , YN , YS ,DY , YL AB , NL Y,NP , ID ) 
DIMENSION  X (80) ,Y(80),XLABI 10) ,YLAB( 10) 

...10-1  GIVES  SINGLE  INCREMENT  SCALE 

''...ID  = 2 GIVES  IW0-1NCREMENT  SCALE 

NX  = XN 

NY  = YN 

XI  = XS 

YI  = YS 

NX  1 = NX  ♦ 1 

NY1  = NY  ♦ 1 

IF(K.NE.O)  GO  TO  30 

DRAW  AND  LABEL  THE  X-AXIS. 

CALL  PLOT ( 0 . , 0. , -3 ) 

DO  10  I = 1 1 NX  1 
CALL  PLOTIO. ,-.07,2) 

CALL  NUMBER (-.07,-.21,.10,X 1,0. ,-l) 

•IF  STATEMENT'  CHANGES  INCREMENT  FOR  ALL  POSITIVE  SCALE. 

GO  TO  (11,9),  ID 
9 I F ( I — i> ) 11,11,12 


11 

XI  = 

XI  ♦ DX 

GO  TO 

i 13 

12 

XI  = 

XI  - DX 

13 

CALL 

PLOT( 0. ,0. ,3  ) 

I F ( I . 

EC .NX  1 ) GO  TO  10 

C/LL 

PLOT ( 1 . ,0 . ,— 2 ) 

10 

CONTINUE 

CALL 

PLOT(— XN,0., —3 ) 

AX  = XN/2 • — . 19* ( 6./7 . ) *NLX/2 • 

CALL  SYMBOL! AX, -.42,. 19, XL AB,0., NLX) 

draw  and  label  the  y-axis. 

CALL  PLOT ( 0. ,0. ,-3) 

DO  20  I = 1 , NY  1 
CALL  PLOT!-. 07,0., 2) 

CALL  NUMBER (-.14, -.07,. 10,YI ,90.,-l > 

YI  = Yl  ♦ DY 
CALL  PL0T(0. ,0. ,3  ) 

1 F ( I .EQ.NY1 ) GO  TO  20 
CALL  PLOT ( 0 • , 1 • ,— 2 ) 

20  CONTINUE 

CALL  PLOT ( 0 . ,— YN , -3 ) 

AY  = YN/2.  — . 14*(6./7. ) *NLY/2 . 

CALL  SYMBOL (-.35, AY, . 14,YLAB, 90 . , NLY ) 

THE  CURVE  MAY  BE  DRAWN  BY  'CALL  LINE*  IN  THE  X-Y  ARRAYS  ARE  FIRST 
FIRST  CONVfRTED  INTO  INCHES.  'CALL  CSCEG'  CONVERTS  AUTOMATICALLY. 
...'ST-30'  SELECTS  THE  TWO-IMCREMENT  OR  THE  ONE-INCREMENT  AXIS. 

30  GO  T0( 31 « 32 ) , ID 

31  CALL  CSEG( 0,K,NP,X, XN»XS,DX,Y, YN , - YS , — OY ) 

GO  TO  33 
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32  CALL  CSEG<  0,K,NP,X  ,%t  , -X S , -DX , Y , YN , -Y S ,-DY ) 

C SET  NPW  uRlOlN  3 ITCHES  BEYOND  PRESENT  PLOT. 

33  CALL  PLUT ( XN  + 3. , 0. ,-3  ) 

RETURN 

t NO 


SUBROUTINE  FbCUT(  1 T , NTH ,NPH I , l X , LNP ) 

L ' SlJu-PbCUT ' PLuTS  THE  PH1=0  AND  PHI=180  CUTS  ON  SAME  PLOT. 

DIMENSION  STT(bu) , SP P ( 80 ) , SF PP C 80 ) 

COMMON  /E VS/PHIN( 20) , ST ( AO) ,S(665) ,SF(665 ) 

1>3  F Ok  M A T ( // 1 X , 1 OH P H 1 , T H F T A = , F 5 • 1 * 1 6 F 7 • 1 / ( 9 X » 1 7 F 7 . 1 ) ) 

U*.  FOKMAT  ( 2X,  SIS) 

I 7u  FORMAT  < Fo. 1 1 3X » I7F7.1/I VX , 17F7. 1 I ) 
no  600  J = 1 , NT  h 
JU  = NTH  ♦ J 
LN  = NTH  ♦ 1 - J 
LL  = LX  - NTH  ♦ J 
ST T ( J ) = -SI  ( LN  ) 

STT(JO)  = ST ( J ) 

SPP(J)  = S ( LN ) 

SPP(JU)  = SUL) 

SFPP(J)  = SF(LN) 

SFPP( JU)  = SF (LL) 
bOO  CONTINUE 

L THF  ARRAYS  ARE  READY  FOR  USE  BY  PBPLOT  RUUTINES  FOR  PLUTT1NC. 

NT2  = 2*N1H 
NT  1 = NTH  + 1 

WklTE(fc,144)  LX.NTH.NPHI ,NT2,NT1 
UP  TO  UlOfbll  It  CMP 

6 1 <J  CALL  PCPLL.T  ( IT  , STT  , 10.,  200.  ,-40.  , 14HTHETA,  DEGRLES,  14, 

1SPP,7. ,2b. ,-4. ,2BHP0WER  IN  OF,  THETA  COMPONENT , 2b , NT2 , 2 ) 
WRITE! 6, 143)  ( STT ( J ) , J= 1 ,MTH ) 

WRI TE(  6,170  PHIN(  I ) , (SPP(J)»J-1 , NTH  ) 

WRI  TE  ( t * l‘,3  ) (SIT  (J)»J=NTl»NT2) 

WRITE (6, 170)  PHIN ( NPHI ) , ( SPP ( J ) ,J=NT1 ,NT2) 

GO  T U oil 

11  CALL  PBPLUT  ( IT,  STT, 10  .,  200.  ,-40.  ,14HThETA,  DEGREES, 14, 

1SFPP ,7. . ,-4. , 26HP0WER  IN  OB,  PHI  COMPONENT ,26»NT2*2) 

WRITE  (c, 143)  ( ST  T ( J ) , J* 1 , NTH) 

WRI TE(t , 17U)  PH  I N( l ), (SFPP(J) ,J=l,NTH) 

WkITt(6,l43)  ( ST  T ( J) , J=NT1,NT2) 

WRITE! b,  170)  pH  IN (NPHI ) , ( SFP P ( J ) , J=NT 1 ,NT2  ) 

1 12  CONTINUE 
RETURN 
END 


BLOCK  OATA 

COMPLEX  ECTT, FCF F, EKTT, ERFF 

Common  /CFLH/EC. TT  ( 66  5 ) , FCFF  ( 665  I / RFL  D/ERTT  ( 66S)  » R ^ F f ( 665  1 
COMPL  TX  ECTT/665* ( I . E- 1 5,0 . ) / , FCFF /665 * ( l.E-15,0. ) / 

C n M p l.  EX  rPTT/665*  (1.  E-15, 0.)/,  ERFF/  66  5*(  l.E-15,0.)/ 

END 
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APPENDIX  C 


EXPANSIONS  FOR  TIIE  COMPUTATION  OF  THE 
BESSEL  FUNCTION  J „ (x) 

In  the  event  that  library  routines  are  not  available  three  different 
expansions  may  be  used  for  the  Bessel  functions  J[(  (x)  for  different  ranges 
of  x and  u . 

For  x<35  and  r>x,  the  power  series  expansion^^ 


Jt.  (x) 


(!) 


GO 

v - (-L2L  ) 

v ! Z-*  k ! ( r + k ) ! '2  1 

k = o 


2k 


(c-n 


is  used.  A more  convenient  form  of  (C - 1) suitable  for  computation  is  given  by 

(JL)" 

2 J 

HB2 


J, 1x1  - n-tn 


(C-2) 


where 


S.  = , 

1 (^  + l)  1 


The  asymptotic  expansion  [ 1 5 ] is  used  to  represent  the  Bessel  function  for 
35<  x < 130  and  » < 3 v/x  . 

1 /2 


\ t \ °° 

J*(x)  = [-tt)  [ cos  (x  - iH1' + r>)  E (- 


1) 


m ( *'  , 2m) 


m = o 


(2x) 


2m 


(C-3) 


- sin 


(*-r«'4')E  ] 

v ' m = o (2x)  J 


where  ( v , m)  is  Hankel's  symbol  defined  by 


m)  = 


, -2m 


m ! 


{<4 


2 2 

1)  (4  yC  - 3*)  . 


.(4  , 2 - (2m-l  )2]  | 


(C-4) 
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The  Meissel  expansion  ^ 1,1 1 valid  for  v real,  large  and  positive 
is  used  to  represent  the  Bessel  function  for  x > 130  and  v >3  </~x.  This 
expansion  is  included  here  for  convenient  reference. 


J ( » + f ) 


1/6 


where 


and 


U 


T (*  ) =1 
o 


UQ ( * ) - ‘ 


T2  (e)  = 


u2(. ) 


77"^ 

+ (? 

it-y) 

X ( v + « 

771 

v (-§-) 

X(  v + t ) 

273 

1 + 
8100 

2 

e 

747 

13 

; + - 

u 


r o 

CD 


E V»  Mr 


r-o 


Tl(<)  = -45 


18 


Uj(.«)  = - 


4 6 

< 

+ 


420  TtT 


283300 


19, 3 _i_5  + ,7 

11340'  810  + 4536 


36 


(C-5) 


Taking  x = v + * and  « = x - v , the  above  expansion  of  the  first  three  terms 
becomes : 


.1 


Jj  (x)  = C|  x 


"C2  X 


/3 


4 50  + — (60  - 1 50  t Z)  + [-2  + 25*2  - 3 5 <4  + 10«5] 


450  - 31  500  x t - 5250  ,2  [ 1 - «2] 


(C-6) 


+ E [52  _ 1900  f2  + 140  t4  - 2500  /’]J 


where 


r (1  / 3 ) 


4 50  * (48) 


176 


/6 


C + (3/4)'  r (2/3) 

2 31  500  "■ 
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